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The fundamental frequencies, anharmonicities, and vibrational modes for the molecules 
HDO, HTO, DTO, and T:0O are calculated on the basis of Dennison’s and Darling’s recent 
analysis of the vibrational spectrum of the HO molecule. The equilibrium constants for nine 
equilibria involving these molecules are given. Experimental values for the reactions HD+H,O 
=H.+HDO and HT+H20=H:2+HTO are compared with those calculated. It is concluded 
that the vibrational potential function is essentially unaltered by isotopic substitution. 





HE recent work of Dennison and Darling’ 

on the vibrational levels of the water mole- 

cule has considerably improved our knowledge 
of the detailed nature of the potential function 
in the ground state.’ As a consequence, an im- 
provement in the calculated levels for the isotopic 
water molecules should result. Several authors?~® 
have given results for HDO and D,O on the older 
analyses of the vibration-rotation spectrum of 
H.O. The agreement between the results of this 
research and those previously published is corre- 


_* This research was completed during tenure of a John 
Simon Guggenheim Fellowship. 

 B. T. Darling and D. M. Dennison, Phys. Rev. 57, 
128 (1940). 

%D. M. Dennison, Rev. Mod. Phys. 12, 17 (1940). 

‘Harald H. Nielsen, Phys. Rev. 62, 422 (1942), has 
recently given values for the frequencies and anharmonici- 
ties slightly different from those of Dennison and Darling. 

he differences are so small, however, as to have little 
effect on the considerations in the present paper. 

* T. Forster, Zeits. f. physik. Chemie B27, 1 (1934). 

*A. Farkas and L. Farkas, Proc. Roy. Soc. (London) 
Al44, 467 (1934). 
(1934) H. Crist and G. A. Dalin, J. Chem. Phys. 2, 735 
(1934) Topley and H. Eyring, J. Chem. Phys. 2, 217 
oar’ Jones and A. Sherman, J. Chem. Phys. 5, 375 

). 


lated directly with the analysis of the H,20O. 
molecule spectrum adopted in most of the papers. 

The discovery of tritium’ (radioactive hydro- 
gen) has made possible further tests of the funda- 
mental assumption of all such calculations, 
namely, that the potential function of a molecule 
shall be independent of the isotopic composition 
of its constituent elements. The new molecules 
which are subject to experimental investigation 
at present are HTO and DTO. The other, T,O, 
has not been obtained in sufficient amounts or 
purity to allow examination, due to the radio- 
activity of tritium (half-life approximately 30 
years’). 


I. FUNDAMENTAL MODES AND FREQUENCIES 


Since substitution of one of the hydrogen 
atoms of H,O by D or T destroys the symmetry 
of the molecule, the form and interactions of the 
normal modes are altered rather seriously (cf. 
Fig. 1). In particular the secular equations for 


( 939) W. Alvarez and R. Cornog, Phys. Rev. 56, 613 
1 ; 

8R. D. O’Neal and M. Goldhaber, Phys. Rev. 58, 574 
(1940). 
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Fic. 1. Normal modes for isotopic water molecules. 


the frequencies must be cubics. Cross and Van 
Vleck® have discussed the problem for the general 
triatomic molecule, arriving at simplified general 
expressions leading to the secular cubic. 
Dennison’s and Darling’s!*" analysis gave as 
the expression for the first-order or simple har- 
monic potential energy of the water molecule 


2 V (ergs) = 8.428 X 105( 67)? 
+8.428 X 10°(67r2)?+0.7676 X 105(Réa)? 
— 2(0.1051) X 10°(6r16r2) (1) 
+2(0.2521) X 105(Réa) (dr167r2)* 


where 6r; and 6r2 are the differences in cm be- 
tween the actual distances between the O atom 
and each of the two hydrogen atoms and the 
corresponding equilibrium distance, R, and 6a is 
the difference in radians between the actual angle 
between the two lines joining the O atoms to 
the hydrogens and its equilibrium value. The 
equilibrium values of R and a are 0.9580A and 
104° 31’, respectively. 

The fundamental frequencies for the six iso- 
topic molecules based on the above potential are 
given in Table I. The hydrogen masses used 
were 1.00813, 2.01473, and 3.01705 for the H, D, 
and T atoms, respectively. 

The normal modes are represented in Table II, 
where 67; always refers to the distance between 
the O atom and the heavier of the two hydrogen 
atoms in the molecule. The dimensionless vari- 
ables, x;, related to the normal coordinates, Q;, by 


xi;= 2x(v;/h)'Q; (2) 


® Cross and Van Vleck, J. Chem. Phys. 1, 357 (1933). 

* Note: The expression (1) differs from that given in 
reference 16 in the sign of the fourth term on the right, 
apparently because of a typographical error in reference 1b. 
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are used. The coefficients are adjusted to give the 
6r1, 6r2, and Réa values in Angstroms. For 
classical vibrations x; ranges from 0 to +Vv2. 

An interesting characteristic of the normal 
modes for the asymmetric molecules is that the 
two stretching vibrations approximate vibrations 
of the hydroxyl groups with respect to the re- 
maining hydrogen atoms, the two frequencies 
being separated because of the different masses. 
This is illustrated in Fig. 1. 

The results of the application of the Teller," 
Redlich" product rule are presented in Table IIT. 
The rule rigorously. requires that insofar as two 
molecules have identical potential energies for 
motions of the constituent atoms they must have 
identical values for the expression 


(xv,)?(axm,)§ 

M*I4IgI¢ 
in which zp; is the product of the frequencies, mm; 
is the product of the masses of the constituent 
atoms, M is the molecular weight, and I4/zJ¢ is 
the product of the principal moments of inertia. 


Column 4 of the table contains these values. The 
check appears to be satisfactory. 


Il. ANHARMONICITIES 


A. Values Uncorrected for Kinetic Energy 
Perturbations 


Dennison and Darling! have shown from ex- 
perimental data that the anharmonicities, X ;;, 


TABLE I. Fundamental frequencies of isotopic 
water molecules. 








. Frequencies (cm!) 





Molecule 3 O1 2 
H.0 3935.59 3825.32 1653.91 
HDO 3883.8 2820.3 1449.4 
DO 2883.79 2758.06 1210.25 
DTO 2830.7 2357.1 1117.9 
T:0 2436.12 2296.63 1017.89 
HTO 3882.6 2365.0 . 1374.5 








in the energy expression 


W= w1(m1+ 4) + w2(m2+ 4) + 03(m3+3) 
+X 11(M1+ 2) (M1 +3) +X 12(m14+ 4) (m2+3) 
+X 20(m2+4) (M2+ 5) +X 23(m2+3) (m3+ 3) 
+X 13(m1 +4) (M3+3) +X 33(a3+ 3) (m3+}) 
10W. R. Angus, C. R. Bailey, J. B. Hale, C. K. Ingold, 
A. H. Leckie, C. G. Raisin, J. W. Thompson, and C. L. 


Wilson, J. Chem. Soc. (London), 971 (1936). 
1Q. Redlich, Zeits. f. physik. Chemie B28, 371 (1935). 
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for the HO and DO molecules are nearly in the 
direct ratio (w;’wj’/wiw;) where the primes refer 
to the H,O frequencies, a result strictly true for 
H,0'* vs. D.O*. Since both the HO and D.O 
molecules are symmetrical, doubts arise that the 
method is applicable strictly to the asymmetrical 
molecules HDO, HTO, and DTO when compared 
with H,O or DO. This point appears to deserve 
further consideration. 

Expressing the potential energy in terms of the 
dimensionless variables x;, 


2V=hyyxPthryexe?+hy3x3?+aix 13+ a2x2' 
+ 3x 3° + 12% 17K 9! + aa1X 2X1 + a13X 7X3 
+ 0431 3°X 1 + A123 2°X 3+ 132% 3°X2+ A123K1K2X3 (4) 
+ Bix1'+ Boxe! + B3x34+ B12% 1x2? + Bigx 1x3" 


+ Bog 22x37 + €1X 1X 29+ €n% 2x13 + 


other odd powers of fourth order, standard 
perturbation theory gives expressions for the 
anharmonicities in terms of the potential con- 
stants. These are given in Appendix I. 

Certain small adjustments in the numerical 
constants in these expressions were made to in- 
sure that the expression for the energy W remain 
a power series in (”;+43). Dennison’s and Dar- 
ling’s' discussion of this point is completely 
applicable here. 

The procedure adopted has been to calculate 
the a; and 6; constants for each isotopic water 


TABLE II. Normal modes of isotopic water molecules. 











ISOTOPIC WATER MOLECULES 


Molecule Modes 

6ri:= 0.0676x, — 0.00647 x. —0.0690x; 

H.O b6r2= 0.0676x, —0.00647x. +0.0690x; 

Ria= 0.000946x, +0.209x» + Ox; 

6r;= 0.0824x, —0.00537x. +0.0063x3 

HDO 6r2= 0.00480x, —0.00559x»2 —0.0965x; 
Ria = —0.00295x; +0.196x2 +0.00104x; 

6ri= 0.0573x; —0.00614x. —0.0591x; 

DO 6r2= 0.0573x1 —0.00614x. +0.0591x; 

Ria= 0.0050x, +0.179x-2 + Ox; 

6ri= 0.0748x, —0.00371x2 +0.0116x; 

DTO bro= 0.01042, —0.00387x. —0.0820x; 
Ria = —0.00895x, +0.172x2 +0.00486x; 

6ri1= 0.0524x, —0.00623x- —0.0542x; 

TO bre = 0.0524x; - 0.00623x.2 +0.0542x; 

Ria= 0.0090x, +0.1638x-2 + Ox; 
6r;= 0.0754x; —0.00479x. +0.00408x; 

HTO 6re= 0.00400x, —0.00509x. —0.0965x3 
Ria = —0.00474x, +0.191 x2 +0.00160x; 
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TABLE ITI. Product rule check on fundamental frequencies. 


=e)" 1 
( M Tal ple 














( mare) 3(v1v2v3)? 
(vivw3)? \ M 











Molecule Talpic 
H.0 a 6.1992 0.7783 
HDO — 2.5204 0.7777 
D,0 a 0.92659 0.7764 
DTO s7isdp «038686 (.7764 
T.0 ona -0:32433.——«0.7776 
HTO ome 1.5930 0.7777 








molecule in terms of those for the H2O molecule 
and substitute them in the above expressions. 
Using the expressions in Table II and equating 
the deformations 67;, 6r2, and Réa for the two 
molecules leads to a linear relation between the 
x; for the isotopic molecule and the x,’ for the 
H.0 molecule, x,’ = oy iijx;. 

Substitution in Eq. (4) gives expressions for 
the a; and 6; for the isotopic molecule in terms 
of those for H2O, a,’ and 8;’. These expressions 
are given in Appendix IT. 

Table IV gives the potential constants, those 
for H.O being those from Dennison and Darling. 

Table V gives the anharmonicities calculated 
from these constants and the equations in Ap- 
pendix I. 


B. Kinetic Energy Perturbations 


These values for the X ;;'s are not quite correct 
because the motions of the atoms contribute zero 
angular momentum with respect to axes defined 
in terms of the equilibrium positions of the atoms 
only at the instant that the atoms are in these 
equilibrium positions. The actual motions con- 
stitute an oscillatory vibrational motion with 
respect to these axes. These energies are to be 
calculated as perturbations of the levels derived 
on the basis of the axes established with respect 
to the equilibrium positions. Weinberg and 
Eckart" have discussed this calculation, and 
Dennison and Darling! describe its application 


122 A. Weinberg and C. Eckart, J. Chem. Phys. 5, 517 
(1937). 



















































































W. F. LIBBY 
TABLE IV. 
Constant H:0 HDO D:0 DTO T:0 HTO 
ay — 322 —272.2 — 197.4 — 203.7 — 149.7 — 208.2 
a2 —47 — 36.9 — 30.85 —22.7 — 24.79 — 34.0 
as 0 431.8 0 264.0 0 433.7 
a12 1 45.66 13.08 14.5 20.58 31.45 
a2} 216 124.7 134.4 101.25 99.1 109.9 
13 0 — 34.39 0 — 67.0 0 —19.5 
a31 — 909 —96.48 — 562.6 — 127.0 — 431.1 —72.4 
23 0 — 126.4 0 —75.9 0 —121.8 
a32 160 74.07 107.9 45.7 86.97 71.8 
123 0 127.2 0 96.74 0 120.9 
By 39 37.91 20.28 25.79 14.02 26.46 
—4.57X 10 e. +0.0123¢. —0.01002e> +1.83X 10~*%es —0.00497e> 
+0.00705«, —0.00553e, +0.0124e, . —0.00824¢e, +0.00944e, 
Bo 2 1.54 1.073 0.908 0.748 1.383 
Bs 35 70.44 18.8 36.43 13.37 71.02 
‘ —2.35X 10 es — 0.00352¢2 —0.00332¢ 
—0.0290e; +0.0375¢e; —0.0535¢e, +0.056¢e, —0.0362e, 
Bie —116 — 72.0 — 61.3 — 46.45 — 42.55 — 56.8 
+0.01002e. —0.0163¢2 0.021 6¢2 —0.0242e. +0.0115¢2 
Biz 212 5.22 112.2 8.62 78.6 4.2 
—0.0212¢. —0.0373¢2 —0.0275¢2 
—0.0140e; —0.0271e; —0.0177e; 
Bos — 108 —97.4 — 57.9 — 54.5 — 40.8 —92.6 
TABLE V. Anharmonicities (uncorrected for KEP). 
H20 HDO D20 DTO T:0 HTO 
Xi — 43.89 —42.27 — 22.65 — 27.84 — 15.71 — 29.43 
+0.0106¢e; —0.0083¢e, +0.01866¢; —0.0124e, +0.0141e, 
Xo — 19.50 4.1* — 15.92 —13.7* — 8.212 —5.91 
X 33 — 46.37 —77.39 — 25.27 — 41.77 — 18.39 —77.70 
—0.00353¢2 —0.00528e. —0.00498¢e2 
—0.0290e; +0.0375e, —0.0535¢, +0.0566¢; —0.0362e; 
X12 — 20.02 —93.7* — 6.76 13.53* —4.27 — 22.45 
+0.01002e> —0.0163¢2 +0.0216¢ —0.0242e2 +0.0115¢€2 
X13 — 155.06 — 16.00 —79.13 —9.01 —56.79 —11.63 
—0.0212¢2 — 0.037 3¢s —0.0275¢2 
—0.0140e, — 0.027 1e; —0.0177e, 
X23 —51.1 — 40.03 — 24.32 — 20.00 — 17.41 — 38.35 
+0.0107e2 +0.0148¢e. 0.0177e2 















to the symmetrical water molecule. The Wilson 
and Howard!* Hamiltonian is used. 

In general for triatomic molecules the first- 
order perturbation energy due to this approxi- 
mate but very simplifying choice of axes is of 
the form (1/2C)P, where C is the moment of 
inertia about the axis perpendicular to the plane 


13E, B. Wilson and J. B. Howard, J. Chem. Phys. 4, 
260 (1936). 


* Values corrected for resonance splitting of 2w2 and w: states. New value for 2w2—w1 is [(2w2—w1)?+a22]4; cf. Section on Resonance Splitting. 


of the molecule and P, is the average square of 
the angular momentum parallel to this axis. 


TABLE VI. Corrections to X;; for kinetic energy 


perturbations (cm). 











H:0 HDO D:0 T20 HTO 
X12 0 4.4 0 0 $5 
X13 0 0 0 0 0 
X23 31.3 14.5 16.2 16.2 11.4 




































































































































ISOTOPIC WATER MOLECULES 
TABLE VII. Anharmonicities (fully corrected). 
4 H.O HDO D:O DTO T:0 HTO 
X11 — 43.89 — 42.27 — 22.65 — 27.84 — 15.71 — 29.43 
—0.00686¢e2 +0.0184e. —0.0150e2 +0.0275e2 —0.00746¢2 
+0.0106¢€; — 0.00836, +0.0186¢; —0.0124e, +0.0141e, 
X22 — 19.50 4.1 — 15.92 — 13.7 — 8.212 —5.91 
X33 — 46.37 —77.39 — 25.27 —41.77 —18.39 , — 77.70 
—0.00353¢2 —0.00528e. —0.00498e2 
—0.0290e, +0.0375¢; —0.0535¢e; +0.056¢, —0.0362e; 
X12 — 20.02 — 89.3 — 6.76 13.53 —4.27 — 18.95 
+0.0100¢. —0.0163e2 +0.021 662 — 0.0242, +0.0115¢. 
2 X13 — 155.06 — 16.00 —79.13 —9.01 — 56.79 — 11.63 
7 —0.0212¢. —0.0373¢2 —0.0275¢2 
—0.0140e, —0.0271e, —0.0177e, 
X93 —19.8 — 25.53 —8.12 —11.3 — 6.21 — 26.95 
‘. +0.0107e2 40.0148. 40.0177 
—0.0081.e, +0.0073e, —0.0155e, +0.011e, — 0.00995, 
DXi; — 76.16 — 61.6 — 39.45 — 21.86 — 25.34 — 42.67 
—0.00272¢2 +0.0005¢€2 —0.0053¢2 +0.0008¢2 —0.0027e2 
TaBLeE VIII. Anharmonicities by product method. conjugate to the Q,’s. Integrating 
H:O HDO DO DTO~ TO HTO * ‘ 
: fv 4 rP ry i,j. id T 
= Xn — 43.89 — 23.86 -—22.81 -16.65 15.80 —16.79 
Xo — 19.50 — 14.98 -1044 — 890 — 7.37 —13.48 
SSM Mey Ee EE Ses gives the remit, Welalhing and Echant® 
: 2 “an an as a =a a |e the result. Weinberg and Eckart' showed 
Xes — 19.81 — 17.12 -1062 -— 961 -— 7.53 —16.23 . 
i3Xij — 761 — 568 —403 —340 —283 —49.15 this to be 
See ii aaa aaa ( ; Vi Ve 
a?(ni +4) (m2+4)| —+— 
TABLE IX. Resonance splitting parameters for (m+2)(m2+3) Pe + » ) 
symmetrical water molecules. as : 
—— h 2 . Vi V3 
(cm) (em-! —_— 2 BS BS ——— 
Molecule (calc) xa. i 2 4 +b (m+ 2) ("s+ 2) ? ) 4 ? 
T V3 V) 
H,0 +71.7 +74.8 
D0 438.5 m 
r,0 +27.4 +e(m+))(m+))( +") 
= \ V3 v2 
TABLE X. Resonance splitting parameters for where a?+8?+c?=1. 
sie nsenenstimeiad motecuies The corrections to the anharmonicities there- 
Molecule (cm™!) calc. fore are h 
‘ 1 . Vy 
: HDO 88.2 en eon Meili aloe 
, DTO 71.8 sania (. ) a ( ¥ ), 
‘ HTO 77.8 c\ee ve 1 
) 
5 7 xy h 2 V1 V3 
, -_ AX 33=— ~) +"), 
Other first-order terms exist but are negligible I\Qe vs 4 
for the water molecules. The second-order terms 
= are negligible also for these molecules. The calcu- AX 1fh\ (™ 4 M3 
. T) . ° 03.=—>— -—— c* —_— — 
, lation of P is made by expressing P., the per- ™ Qc\ 29 - a , 
5 pendicular component of the “internal angular 
4 momentum,” as (—aQ2—cQ3)Pi+(aQi+6Q3)P2_ if the constant term is neglected to preserve the 


+(cQ:—bQ3)P3, where the P;’s are the momenta form of the expression for the energy as a power 
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Fic. 2. Theoretical values for (HDO)(H2)/(H20)(HD). 


TABLE XI. Frequencies and anharmonicities of 
hydrogen molecules. 









































He HD D2 DT Te HT 
4405.30 3817.09 3117.07 2845.64 2546.50 3598.14 
62.66 47.01 31.36 26.14 20.94 41.80 
TABLE XII. Ho° (calories per mole). 
Molecule H20 HDO D:0 DTO T:0 HTO 
Ho® 13,241 11,480 9,682 8,951 8,148 10,771 
Molecule He HD De DT Te HT 
H,° 6,253 5,423 4,433 4,049 3,625 5,114 
TABLE XIII. (F—Hp°)/T values (calories/degree). 
Molecules 25°C 100°C 200°C 300°C 400°C 500°C 
H.O —37.20 -—39.00  -—40.88 —42.40 —43.70 —44.90 
HDO —39.80 -—41.55 —43.50 —45.10 -—46.40 —47.50 
D:O —39.42 —41.15 -—43.10 -—44.65 -—46.00 —47.20 
DTO —41.55 -—43.20  -—45.20 -—46.80 -—48.15  —49.40 
TO —41.02 —42.64 -—44.58  -—46.20  —47.57 —48.72 
HTO —40.45 -—42.30 —44.10  -—45.70 —47.05 —48.30 
He —24.37 -—25.90 —27.48 -—28.78 -—29.85  —30.88 
HD —27.48 -—29.06  -—30.60 -—31.88 -—33.00 -—34.00 
De —27.72 —29.26 —30.80 —32.15 —33.25  —34.25 
DT —30.07 31.64 —33.35 —34.60 -—35.67 —36.65 
Te —29.70 —31.26 —32.82 -—34.18  —35.26  —36.20 
HT —28.55 -—30.10 -—31.66 —32.94 -—34.08 —35.03 








series in (n;+3). This point was discussed in a 
preceding paragraph. 

Table VI gives the corrections and Table VII 
the final values for the X ;;. 

Table VIII gives the X;; values calculated by 
the method used by Darling and Dennison! for 
D.O. The agreement for individual X;;’s is quite 
good for the symmetrical molecules and poor for 
the asymmetric ones, though the contributions to 
the zero-point energy, }>.X,;, are in fair agree- 
ment in all cases. The coefficients €; and ¢€2 are 
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Fic. 3. Theoretical values for (HTO)(H2)/(H2O)(HT). 


taken to be of the same order as the other quartic 
coefficients, the 6’s. With this assumption their 
effects on }>/X,; are of the order of 1 cm" 
which is quite small. More accurate experimental 
investigation will require their evaluation, of 
course. Since the e’s refer to bending coordinates, 
considerations of the type used by Redlich'™ 
further justify the expectation that they are 
small. 


III. INFRA-RED SPECTRA AND RESONANCE 
INTERACTIONS 


A valuable check on the results discussed 
above would be afforded by an experimental de- 
termination of the infra-red and Raman spectra 
of HDO by subtraction of the HxO and D,O 
spectra from the composite obtained for a 50 
percent solution of H,O and D.O. 

Fuson, Randall, and Dennison'® have given 
values for about eighty lines of the infra-red 
spectrum of HDO obtained in this way. An 
attempt was made to check the frequencies and 
anharmonicities calculated above against these 
lines. All lines definitely known to belong to 
HDO were attributable to relatively probable 
transitions in the calculated infra-red spectrum. 
The fit usually was as close as 1 or 2 cm. In the 
absence of intensity data, however, a completely 
conclusive test can hardly be obtained with so 
few lines. 

The expression (3) for the energy levels is 


144Q,. Redlich, J. Chem. Phys. 9, 298 (1941). 
15 N. Fuson, H. M. Randall, and D. M. Dennison, Phys. 
Rev. 56, 982 (1939). 
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directly applicable, except for the Fermi!® reso- 
nance effect arising in the cases of certain over- 
tones. Dennison and Darling have discussed this 
problem for the symmetrical HO, and D.O 
molecules. They found serious interaction be- 
tween the levels (73, 21, m2) and (m3+2, ,—2, n2). 
The levels are split symmetrically, the new sepa- 
ration being (67+ 6,7)! where 69 is the separation 
for no resonance interaction and 6; is given by 


N3,N1,N2 


61 => 2bns +2, 21-2, n2 


=hc[n1(m1—1)(m3+1)(n3+2) }by 
where 

2 
Qi3  &1@31 1232 12032 


ee i mas Wieted 


WI 2w1 
TABLE XIV. Masses and moments of inertia. 
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Table IX contains y values for the sym- 
metrical molecules. 

For the molecules HDO, HTO, and DTO the 
most seriously resonating levels are (m3, 11, M2) 
and (3, m1—1, m2+2). For these 


5; =he[my(n2+1)(m2+2) }by 
where y= a2;/V2. 
Table X gives these y values. 
IV. ISOTOPIC EQUILIBRIA BETWEEN THE WATER 
AND HYDROGEN MOLECULES 


The equilibrium constants for the following 
equilibria can now be evaluated. 


(1) HOO+HD=HDO+H, 
(2) HXO+HT=HTO+H, 











Molecule Mass 1 ABc *10'29 (g cm?) (3) D,O+ DT = DTO+ D, 
H;O 18.01626 5.869 (4) DOO+HD=HDO+D, 
DO 20.02946 40.67 = 
T.0 22.03410 120.5 (5) H:0+D:0=2HDO 
HDO 19.02286 16.6 (6) T,0+DT=DTO+T, 
DTO 21.03178 70.84 = - 
HTO 20.02518 29.40 (7) T:O+HT=HTO+T: 
I 0.104° (8) H.0+T.,0=2HTO 
H. 2.01626 0.4599 9 > -—_ 1 
Ds 4.02946 0.9191 (9) D:O+T:0=2DTO 
T 6.03410 1.376 . — 
HD 3.02286 0.6131 The fundamental frequencies and anharmonici- 
DT 5.03178 1.102 a , . . . ‘ven i 
HT 4.00518 0.6804 ties for the hydrogen molecules are given in 
Table XI. 
TABLE XV. Equilibrium constants. 
Equilibrium 20°C 100°C 200°C 300°C 400°C 500°C 
(HDO)(H2) 
1) ——-— 3. 2.65 2.09 1.75 1.54 1.37 
() G.0)(HD) " 
(HTO)(H2) 
—__——— ’ 3.83 2.66 2.09 Ly 1.52 
(2) (H,0)(HT) 6.24 
(DTO)(Dz) 
So ar 1.37 Laz 1.15 1.12 1.08 
(3) (,0)(DT) 1.58 
(HDO) (D2) 
—_ e 1/2. 1/1.78 1/1.55 1/1.40 1/1.29 
© aidID) 1/3.00 /2.27 / / / / 
(HDO)? 
—— 4.10 4.13 4.16 4.17 4.17 4.1 
©) .0)(D0) : 
ie, (DTO)(T2) 
(6) (T:0)(DT) 1/1.77 1/1.45 1/1.36 1/1.24 1/1.19 1/1.18 
(HTO)(T2) 
(7) (T.0) (HT) 1/5.00 1/3.10 1/2.27 1/1.84 1/1.62 1/1.52 
(HTO)? 
(8) (H.0)(T.0) 18 3.57 3.71 3.77 3.82 3.83 
(DTO)? . 
(9) (D.0)(T0) bY 3.66 3.74 3.78 3.78 3.79 








16 FE. Fermi, Zeits. f. Physik 71, 250 (1931). 

















TABLE XVI. 
Temp. (°C) Krfheor Kigxp Reference 
20 3.78 2.67 (1) 
5.2 (2a, b) 
Corrected 
from value 
of 3.8 for 
liquid water 
527 137 1.28 








1 A. Farkas and L. Farkas, Trans. Faraday Soc. oe 1076 (1934). 
2a Bonhoeffer and Rummel, Naturwiss. 22, 45 (1934). 

2K. F. Bonhoeffer, Zeits. f. Elektrochemie 40, 469 (1934). 

3 See reference 4 in text. 


The values for D2, and HD are those given by 
Teal and MacWood.!? The other values were 
calculated from the He values by use of the 
ratios of the reduced masses of the molecules. 

Table XII contains the H,° values for the water 
and hydrogen molecules calculated from the 
fundamental frequencies and anharmonicities 
contained in Tables I, VII, and XI. 

Table XIII contains (F—H,°)/T values for the 
water and hydrogen molecules, nuclear spin 
entropies being omitted. 

The masses and moments of inertia used are 
listed in Table XIV. 

The equilibrium constants are given in Table 
XV. K, and K» are plotted in Figs. 2 and 3. 
They should be accurate to about two percent. 


V. DISCUSSION 


The experimental values for (HDO)(H2)/ 
(H20)(HD) reported in the literature are given 
in Table XVI. 


7G K. Teal and G. E. MacWood, J. Chem. Phys. 3, 
760 (1935). 


TaBLe XVII. Experimental values for 
(HTO)(H2)/(H20)(HT). 











Temp. (°C) K Theor. KExp. Source 
20 6.24 6.3+0.08 Taylor and Black 
Private communi- 
cation 
28 5.80 5.3+0.5 Author 
280 2.00 1.9+0.1 Author 








The analogous tritium constant (HTO)(He2)/ 
(H.O)(HT) can be measured by equilibration 
with Pt of a H2O, He system containing tritium 
in both the H2O and He species as HTO or HT. 
The concentration of tritium will be so low 
(~10-"M) that no precautions need be taken 
about other equilibria. After equilibration the 
hydrogen and water vapor are introduced into 
the gas of a Geiger-Miiller counter or electroscope 
and the relative specific radioactivities deter- 
mined. The author has made a few measurements 
of this constant and a more extensive study is 
being made by Mr. James F. Black and Professor 
H.S. Taylor of Princeton University. Table XVII 
gives the results to date. 

In view of the agreement for the tritium case 
and the difficulties of the determination in the 
deuterium case due to other equilibria, the author 
believes that the fundamental assumption of the 
essential identity of the potential functions of 
isotopic molecules is correct. 

The author is much indebted to Professor H. 
S. Taylor and Mr. James F. Black for permission 
to quote their data. He wishes to thank the John 
Simon Guggenheim Memorial Foundation also 
for the Fellowship during which this work was 
completed. 


APPENDIX I. EXPRESSIONS FOR THE ANHARMONICITIES IN TERMS OF 
THE POTENTIAL CONSTANTS 
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APPENDIX II. EXPRESSIONS FOR THE POTENTIAL CONSTANTS OF ONE 
ISOTOPIC MOLECULE IN TERMS OF THOSE OF ANOTHER 
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The Surface Recombination of H Atoms and OH Radicals 


W. V. SMITH 
Jefferson Physical Laboratory, Harvard University, Cambridge, Massachusetts* 


(Received November 20, 1942) 


A diffusion method of investigating surface recombination of atoms by substances of low 
catalytic power is presented. By this method we have evaluated the coefficient of recombination 
y for H and OH on Pyrex for a temperature range of over 500°C. For smaller ranges of temper- 
ature we have found the recombination of H on a number of salts, and have found this recom- 
bination to be strongly dependent on the dryness of the surface. All salts investigated save 
KCI strongly recombined H when dry. Salts recombining H could not be investigated for 
OH recombination. However KCI, which did not recombine H, did recombine OH strongly. 
A number of strong dehydrogenating and dehydrating agents were tested for differences in 
H and OH recombination, but no difference was found. 





INTRODUCTION 


LTHOUGH the investigation of surface 
recombinations of many atoms and radicals 
is still in the qualitative stage, quantitative 
results have already been obtained for the rate 
of recombination of H atoms on glass! and on 
certain metals.? By a modification in the experi- 
mental arrangement of previous investigators in- 
volving a change in the theory of the experiment, 
we have been able to extend these data and 
especially to find the variation of the coefficient 
of recombination y (the fraction of times an atom 
striking the surface leaves as part of a molecule) 
with temperature. 


1. EXPERIMENTAL PROCEDURE 


The experimental procedure consists in allow- 
' ing atoms that are formed by a low pressure 
gaseous discharge to diffuse through a side arm 
part of whose inner surface is coated with the 
material under investigation. Some of the atoms 
are lost to the coated walls and the remainder 
which succeed in passing through the side arm 
are detected by the heat of recombination they 
liberate at a small probe coated with some highly 
catalytic substance (Pt to detect H?; KCI to 


* Now at Massachusetts Institute of Technology, Radia- 
tion Laboratory. 

1H. G. Poole, Proc. Roy. Soc. 163, 404 (1937); W. 
Steiner, Trans. Faraday Soc. 31, 962 (1935). These values 
were obtained by complicated computations under con- 
ditions where the surface recombination was only a correc- 
tion to a predominant volume recombination. 

2K. F. Bonhoeffer, Zeits. f. physik. Chemie 113, 199 
(1924). 
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detect OH).’ A high probe temperature indicates 
that few atoms are lost to the walls of the side 
arm, whereas a low probe temperature shows 
that many atomsare lost to the walls. By varying 
the probe distance, the decay in atomic concen- 
tration with distance from the juncture of dis- 
charge tube and side arm can be found, yielding 
an absolute measure of the coefficient of recom- 
bination; and with probe distance fixed, vari- 
ation of the temperature of the coated section of 
the side arm yields the temperature dependence 
of the recombination coefficient. 

Figure 1 is a diagram of the apparatus. The 
hydrogen was evolved electrolytically from a 
concentrated solution S of NaOH and passed 
through a water trap T which regulated the 











ol___»___p 
| 


region ~ region x= 











Fic. 1. Diagram of apparatus. S and W are sources of 
H2 and H.O, respectively. T and 7” are traps. C and C’ are 
capillary tubes restricting gas flow through discharge 
tube D. A is auxiliary discharge. P is probe at distance X 
from discharge. The temperature of region 1 of length L is 
varied. Region 2 is at room temperature. 


3H. S. Taylor and G. I. Lavin, J. Am. Chem. Soc. 52, 
1910 (1930). 
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percent of water vapor in the gas.‘ It was thence 
pumped through a capillary tube C and on 
through the discharge tube D by a Cenco Hyvac 
pump. A continuous flow of hydrogen through 
the discharge tube ensured freedom from im- 
purities. The water vapor source was simply a 
bulb W of distilled water, the vapor of which 
was pumped through a capillary tube C’ into 
the discharge tube. The trap 7”, placed in dry 
ice, kept the water vapor out of the pump. A 100- 
milliampere current passed through the tube dis- 
sociated the gases respectively into H+H and 
H+OH*>. 

The gases were investigated at pressures less 
than 0.1 mm, the pressure being governed by the 
capillary restriction and the pumping speed. The 
dark space of an auxiliary discharge A with 
plane electrodes was used to measure the pres- 
sure. It was calibrated against a McLeod gauge 
for hydrogen and against the vapor pressure of 
ice at known temperatures for water vapor. All 
stopcocks were lubricated with special low vapor 
pressure grease and were placed out of range of 
the atomic hydrogen. 

The temperature of the probe P was measured 
by a copper-advance thermocouple with the 
cold junction 2 or 3 cm from the hot junction, 
thereby measuring the temperature difference of 
the gas and the heated probe. The whole thermo- 
couple was separated from the low pressure gas 
by a glass wall, thin at the tip, so that the 
thermocouple junction (at atmospheric pressure) 
assumed the temperature of the platinum 
coating on the low pressure side of the glass. 
The probe was glass-blown to the side arm to 
avoid waxes or greases which might be attacked 
by atomic hydrogen. 

The experimental procedure is based on the 
following considerations: The surface recom- 
bination of atoms must compete with recom- 
bination by triple collisions in the gas phase. 


4 Addition of about two percent water vapor poisons the 
glass walls to H recombination. [R. W. Wood, Proc. Roy. 
Soc. 97, 455 (1920).] Other poisoning agents such as 
metaphosphoric acid were avoided as they contain an 
indeterminate amount of water vapor and also make 
cleaning of the apparatus more difficult. 

5 Spectroscopic observations have shown H and OH to 
be the only important decomposition products of an 
electric discharge through water vapor. [K. F. Bonhoeffer 
and T. G. Pearson, Zeits. f. physik. Chemie. B-14, 1 (1931).] 








Simple calculations show that with y of the order 
of 10-5 (the smallest value found was 210-5), 
in apparatus of the dimensions used a pressure 
of 0.1 mm is sufficiently low to make volume 
recombination negligible. However at this low 
pressure diffusion becomes very rapid. Therefore, 
one cannot very well observe the recombination 
by the decay of concentration along a tube 
through which gas is flowing. A pumping speed 
high enough to make diffusion negligible is not 
practicable, particularly for H. Furthermore, the 
smaller y, the smaller is the pumping speed 
required to produce a given concentration 
gradient of atoms along the side arm. Thus, 
rather than attempt cumbersome and _ inac- 
curate. corrections for the influence of diffusion 
on flow velocity, it is far simpler to rely on dif- 
fusion alone to transport the atoms to the probe. 
It is in this reliance on diffusion rather than 
pumping that our experiment differs from pre- 
vious ones. 

Adoption of a diffusion scheme however 
requires that our detecting apparatus influence 
the atomic concentration as little as possible. 
Two techniques are available which have less 
disturbing influence than the one we adopted. 
They are the spectroscopic method of Steiner’ 
and the Wrede slit method.* Both techniques are 
rather complicated however, the Wrede one 
especially so, and both have difficulties in inter- 
pretation, particularly the Steiner spectroscopic 
method. Simplicity of design was a consideration 
in our experiment as we investigated numerous 
different samples, necessitating frequent cleaning 
of our apparatus. Furthermore, we are able to 
show that most of our results are unaffected by 
the recombination at the probe. 


6 Steiner (reference 1), reports volume recombination 3 
times as effective as surface recombination for a pressure 
of 0.3 mm and tube diameter of 25 mm. This is not incon- 
sistent with our criterion when due account is taken of 
his different pressure and diameter. Furthermore, his lower 
value of y indicates better poisoned glass than that with 
which we worked. 

7Z. Bay and W. Steiner, Zeits. f. physik. Chemie B3, 
149 (1929). The relative atomic concentration at points 
along a side arm from the discharge tube is deduced from 
the spectra of weak localized high frequency discharges 
induced at these respective points. 

8 This method relies on the different velocities of effusion 
of molecules and atoms through a small slit. The differ- 
ences in theoretical slit sizes necessary and the sizes used 
by experimenters has been criticized by Poole [Proc. Roy. 
Soc. 163, 424 (1937) ]. 
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Besides simplifying the mathematics, the dif- 
fusion method avoids stray discharges as there 
are no grounded points along the side arm. We 
wish to emphasize this point as it is particularly 
important if any significance is to be attached 
to low probe temperatures, which might other- 
wise be ascribed to the effect of stray discharges. 
With a flow through the side arm, which always 
ends with the pump as a grounded point, it is 
difficult to be sure stray discharges are avoided, 
especially when an alternating voltage is applied 
to the discharge. 


2. THEORY OF THE EXPERIMENT 


The atomic concentration at various distances 
from the discharge along a side arm will vary, as 
different parameters such as pressure, tem- 
perature, and degree of dissociation in the 
discharge are varied. The nature of this variation 
will differ according to the mechanism responsible 
for the recombination process. Four types of 
mechanisms are possible. They are: 

1. A first-order wall effect with total recom- 
bination directly proportional to atomic con- 
centration. Such an effect results from the 
assumption that each atom striking the wall has 
a definite probability of sticking to the wall and, 
once it sticks, is certain eventually to combine 
with a second atom and leave the wall as part 
of a molecule. This might be the mechanism of 
recombination on poor catalysts. The assumption 
of the wall being completely covered with a 
monolayer of atoms, independent of pressure, 
also results in first order recombination as in this 
case the number of molecules leaving the wall 
per second can vary only in direct proportion to 
the number of atoms striking the wall per second. 
Good catalysts, such as Pt, probably fall in this 
class. Whatever the mechanism responsible for 
such an effect, it can be described by assigning a 
coefficient of recombination y to each atom 
striking the wall signifying the probability that 
the atom will eventually leave the wall as part 
of a molecule. 

2. A second-order wall effect. Recombination 
proportional to the square of the pressure results 
from the assumption that each atom striking the 
wall adheres to it a short time only and that 
recombination will occur if, during this time, 
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another atom strikes the wall within acertain 
(small) distance of the first atom. Any practical 
case intermediate in pressure dependence between 
1 and 2 can be resolved into these two parts. 

3. Triple collisions in the gas phase. 

4. Preponderant effect of the probe. If there 
were no wall recombination nor any losses of 
atoms by triple collisions, but only a large recom- 
bination at the probe, there would still be a 
variation of probe temperature with distance 
from the discharge as the transport of atoms to 
the probe is a function of the distance and the 
coefficient of diffusion. We cannot deduce from 
the small probe area (5X10-* cm?) that it is 
negligible, as the recombination is 10* or more 
times as great on Pt as on glass. (See Table IT.) 

From our measurements we want to decide 
between cases 1, 2, 3, and 4. For this purpose, 
we carry through the theory for each of the four 
cases and find which checks best with experiment. 


Let »=number of atoms/cc 
x=distance along side arm from junction 
with discharge tube 
D=coefficient of diffusion 
F=drift velocity of atoms along side arm 
y=coefficient of recombination 
S=cross-sectional area of side arm 
C=circumference of side arm 
d=average atomic speed 
p= pressure. 


The subscript 0 in the following derivations 
indicates the junction of the side arm with the 
discharge tube (Fig. 1). 

We will assume the concentration is constant 
over the cross section so that n=n(x) only. We 
assume also F= F(x).° Hence for all four cases 


* For the equation of diffusion to be applicable to our 
problem, we must have 2r>>A. (r is the radius of the side 
arm; A, the mean free path.) Also, will equal n(x) only 
if the diffusion is rapid enough to wipe out concentration 
gradients within any cross section which the losses to the 
walls tend to establish. A typical time for diffusion to 
wipe out these gradients is given by the relaxation time 
for diffusive equilibrium in an enclosure of dimensions r: 
t=r?/D=3r?/dd. This time must be small compared to the 
time necessary for destruction of all atoms in the volume V, 
neglecting the effect of diffusion. It is 


t<—” ae 
~ n0Sy/4 07 





Therefore our complete condition for validity of our 
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the equation of diffusion applies in its simplest 
form: 
dn 


nF = —D—. (1) 
dx 


A second equation between n, F, and x is 
obtained by expressing the decay of the concen- 
tration along the length x of the side tube in 
terms of the mechanism by which the atoms are 
removed. We may derive the dependence of n 
on x for cases 1 to 3 simultaneously, assuming 
constant temperature along the side arm. We 
first equate the number of atoms lost to the wall 
per cm per second to the gradient of atomic flow 
through the tube [—S(d/dx)(nF)]. The losses 
to the walls for the three cases are: 


: for first-order wall re- 
(a) —97C combination, 
(b) aCn? for second-order wall (2) 


recombination, 


(c) bpSn*+cSn* for volume recombina- 
tion. 


In (c) the first term is for molecules as third 
bodies and the second for atoms. As our partial 
pressure of atoms is in general small compared 
to the total pressure, we will ignore the second 
term.!° We then have 


wale 

ae 

d 4 

— S—(nF) =: (3) 
dx aCn? 


bpSn?. 





Combining this equation with (1) and solving 
subject to the boundary condition that (dn/dx) 
=0 at n=0 [from Eq. (1), dn/dx is proportional 
to m and in particular approaches 0 as m ap- 
proaches 0 since although F is variable, it remains 


assumptions is \/2<r<2X/3y. In practice we have \=0.2 
cm, r=0.55 cm, and y=107' to 10-7. We see therefore 
that our assumptions »=n(x) and F= F(x) are valid. 

10 The main conclusion of this derivation—that n/no 
is a function of mo, approaching 1 for small mo, would 
still hold if this second term were included, but eliminating 
it simplifies the mathematics. (Data in the literature are 
not in agreement on the relative effectiveness of atoms 
and molecules as third bodies promoting volume recom- 
bination.) 


finite ] we have 
(A exp (6Cy/4SD)'x+B exp —(iCy/4SD) x 
No 


n= [1+(naC/6SD)'x }? (4) 





No 


[1+(nobp/6D)ix 2 








In the expression for first-order wall recom- 
bination, with a tube of finite length both A and 
B+0. The negative exponent represents a com- 
ponent decreasing with increasing x. The positive 
exponent represents a stream of atoms reflected 
from the end of the tube and therefore increases 
with increasing x. We may however use the rule 
of thumb that for values of ”/mo less than about 
i (a condition satisfied by our experimental 
results) the atoms reflected from the tube end 
may be ignored; hence, A1 and B=0. With 
obvious changes of constants, we therefore have 


(a) n=mno exp — (By'x) 
for first-order wall recombination 
with B=(6C/4SD)}, 
No (5) 


n= 
[1+ XK (or K’)no}x }* 





(b) 


for second-order wall recombination 


(or for volume recombination). 


We note that 8 is proportional to D-} and hence 
to pi as the coefficient of diffusion D is inversely 
proportional to the pressure. 

To estimate the maximum influence of the 
probe alone on the concentration (case 4), we 
will make the assumption that all the atoms 
striking the cross section where the probe is 
located are eventually recombined at the probe. 
We will first analyze the problem under condi- 
tions of constant temperature along the side arm. 

At X we have destruction of atoms at a rate 
(nd/4)S. These atoms are being supplied by dif- 
fusion down the side arm at a rate SnF. There- 


fore at X 
F=3/4. (6) 


In the region from x=0 to x=X, where there 
are assumed to be no losses of atoms by recom- 
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Fic. 2. Variation in the ratio of temperatures of two 
probes plotted against the temperature of one of these 
probes. Discharge in He; pressure=78u. 1= Theoretical 
curve for loss of H by three-body collisions. 2 = Theoretical 
curve for first-order wall recombination of H. 3 = Observed 
data. (The constants of the theoretical curves were chosen 
to fit the observed data at the point of intersection.) 


bination, we must have conservation of mass 
transport along the tube, hence 


d 
—(nF)=0. 
x 


(7) 


Combining (7) with the equation of diffusion 
(1) we get 


(8) 


This means that there is a linear fall of atomic 
concentration along the side arm furnishing a 
gradient to transport atoms to the probe by 
diffusion. From Eqs. (1), (6), and (8), we there- 
fore derive. 


n=ny—Kx/D, where K is some constant. 


K @6 n(X)d 
(F)x= =-, whence K= ” 
n(X) 4 
Substituting this value of K in (8) and evalu- 
ating at x=X, we have 


X= where G=—. (10) 

4D 1+GX 4D 
Assuming that not all the atoms crossing X 
recombine on the probe, but that a large fraction 


do so (the condition in actual practice) has the 
effect of making G smaller. We will therefore 





(9) 








n(X)=no— 
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take the magnitude of G as unknown and examine 
only the dependence of our equation on X, 
pressure, and temperature. The latter will prove 
most useful for correlation with experiment. For 
this purpose, we will want a solution for the case 
where the temperature of a portion of the side 
arm far from the probe is varied (thus changing 
the density and the coefficient of diffusion). The 
calculation for this case is straightforward though 
a bit tedious. The resulting formula is 


No 


| i ae 
1+6(x-14+——1) 
TD 


Here 7; and D, are temperature and coef- 
ficient of diffusion in region 1 and J» and Dp are 
the corresponding values in the discharge tube 
and in region 2. (See Fig. 1.) 

To interpret our results, we must relate probe 
temperature to atomic concentration. Let us 
define ‘‘probe temperature” t= 7,—T where 7, 
is the actual equilibrium temperature of the 
probe and 7» is the temperature of the side arm 
in which the probe is located. We assume 


(11) 





n= 


nv 
—Agq=hkt, (12) 
4 

where A is the probe area, g the heat generated 
per atom striking the probe, and k the conduc- 
tion-convection constant. Heat loss of the probe 
by radiation [which would introduce a term 
k'(T,—T»)*] is neglected. 

A complication in our investigation is the 
necessity of assuming this linear relation between 
probe temperature and atomic concentration, as 
we cannot divorce the investigation of loss of 
atoms along the tube from a simultaneous inves- 
tigation of probe response. Such an assumption 
is reasonable however, particularly in the case 
of recombination of H on Pt as it is known that 
metals are quickly covered with a monolayer of 
atoms at temperatures ranging from below room 
temperature to several hundred degree centi- 
grade. If this monolayer is maintained over the 
partial atomic pressure range investigated, the 
number of molecules leaving the surface per 
second must be proportional to the number of 
atoms striking the surface per second. In any 
case, a small deviation from linearity in response 
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would not be serious as most of the concentration 
variations are exponential in character. 


3. RECOMBINATION OF H ON PYREX— 
DISCHARGE IN H: 


All the samples of Pyrex used were cleaned by 
scrubbing and rinsing in distilled water, as pre- 
liminary investigations showed no difference in 
results even at high temperatures using other 
methods of cleaning such as washing in KOH, 
KHCrO,, H.SO,, or HF. 


A. Check on Three-Body Collisions and 
Second-Order Wall Effects 


In order to distinguish between Eqs. (5a) and 
(Sb), we note that (5a) predicts 2/m» independent 
of mo, whereas (5b) predicts /no strongly de- 
pendent on mp, and indeed approaching the 
value 1 as mo approaches 0. The following experi- 
ment was therefore performed. Two parallel side 
arms of equal diameter were attached to the 
discharge tube, terminating in probes respec- 
tively 37 cm and 150 cm from the discharge. The 
atomic concentration was varied over wide 
ranges without materially changing the pressure 
or the coefficient of diffusion by simply varying 
the current through the discharge tube. The 
ratio of the two probe readings was then plotted 
against the temperature of the probe closest to 
the discharge (Fig. 2). The observed approximate 
constancy of this ratio, at a value far from 1, is 
ample evidence that /mo is not dependent on mo. 
We conclude therefore that neither three-body 
collisions nor second-order wall recombination 
on moist Pyrex walls are important factors in H 
recombination at the pressure and temperature 
investigated (78 microns, and 23°C). 


B. Check on the Effect of the Probe 


If the probe recombination predominates, a 
change of temperature of region 1, Fig. 1, should 
have an effect predicted by Eq. (11). This is 
checked in Fig. 3 where the probe temperature, 
indicating the number of H atoms arriving per 
second, is plotted against the temperature of 
section Z of the side arm. For comparison the 
theoretical curves derived from Eq. (11) (pre- 
dominant probe recombination) are also plotted. 
These theoretical curves are adjusted to fit the 
experimental values at room temperature. 





Obviously at temperatures other than room 
temperature the probe recombination does not 
predominate. However as the surface recom- 
bination is lowest at room temperature (see 
Fig. 6), this argument does not eliminate the 
possibility of important disturbance from the 
probe at room temperature. Values of y calcu- 
lated by ignoring the probe effect may however 
be considered maximum values, and all results 
for various surfaces at various temperatures 
yielding a y considerably larger than the least 
of these values may be considered reliable. Fur- 
thermore, the assumption of negligible probe 
effect appears justified even at room temperature 
in view of the agreement between experiment and 
theory shown below in Figs. 4 and 5. This result 
is plausible since the pressure is kept so low and 
the probe area so small that we should expect the 
first-order wall effect to predominate. 


C. Computation of y at Room Temperature 


1. Computation from the Decay in Atomic H 
Concentration Along the Side Arm. (Fig. 4) 


From Eggs. (5a) and (12) we have 
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t= exp — (By'x). (13) 
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Fic. 3. Variation of probe temperature with temperature 
of furnace located at region 1 of Fig. 1. H»2 pressure = 80u. 
1 = Theoretical curve for atomic gradient caused by recom- 
bination at probe, normalized to observed data at 20°C. 
2=Observed data. A: L=58; X =70. B: L=20; X =50. 
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Fic, 4. Variation of probe temperature with distance 
from a He discharge. Pressure =75y. Circled points corre- 
spond to circled point of Fig. 5. 




















Therefore in a series of experiments in which mo, 
d, etc., are kept constant, we derive 


log. t= —Boy'x+C’, 


where C’=some constant. 

For p=75yu, Bo=3.73 (see Section D). 

Three lines are drawn through the points in 
Fig. 4 showing the mean and extreme values of 
the slope. The corresponding values of y cal- 
culated are 


(14) 


0.75 
y=2.01- 10-5. 
4.56 


The variations in y are of the order of 2 or 3 
from the mean. As numerous different samples of 
Pyrex were used the variation in y is reasonable. 


(15) 


2. Computation from Pressure Variation of no/n 


Equation (5a) may be written 


B’xoy' 
ai } 
log. (mo/n) ( _ )e 


where # is in microns and 8’ is the value of 8 at 
75 pressure. 

In Fig. 5, log. 2(30)/n(145) is plotted against 
p'. The theoretical curve 16, a straight line 
through the origin, gives a good fit to the experi- 
mental points. Some deviation is observed at 
pressures below 80 microns. This is possibly 
caused by a breakdown in the assumption of 
uniform drift velocity F along the tube cross 
section since for these pressures the mean free 





(16) 
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path (1.1 mm at 80yz) is no longer negligible com- 
pared to the tube diameter (1.1 cm). 
The value of y calculated from 16 and Fig. 5 is 


y=1.05X10-5, (17) 


The probe temperatures at 30 and 145 cm used 
in Fig. 5 for a pressure of 75y are also plotted on 
Fig. 4 (the circled points). They are seen to 
represent a value of y nearer the lowest value 
than the mean value deduced from Fig. 4. These 
values (see 15) are 0.75X10-5 and 2.01*10-5. 
The check of 1.05 with 0.75 is exceilent. 


3. Comparison with Other Values 


Steiner! reported y on moist glass varying 
between 2X10-* and 9X10-*. Poole! reported a 
value of y of 1.6X10-*. Our values lie between 
these two extremes. Poole’s results are for glass 
under the actual conditions of the discharge 
however, for which y is to be expected com- 
paratively large. The agreement of these values 
with ours is very good in view of the complicated 
mathematics used in the work of the other two 
investigators where the wall recombination was 
a minor effect. The simplicity of interpretation of 
our results plus the fact that the wall recombina- 
tion was the preponderant agent removing the 
atoms in the pressure region we investigated 
would seem to give greater confidence to our 
results. 


D. Experiments at Other Temperatures 


In these experiments, the pressure, tube 
current, and distance X are kept constant while 
the temperature of a section L near the discharge 
is changed by means of a furnace or low tem- 
perature bath (see Fig. 1). In some runs a parallel 
check side arm and probe was used to see if the 
atomic H production remained constant during 
the run. It was found that the probe temperature 
in this check side arm would vary by a factor of 
2 or 4 for the first few hours and then remain 
constant often to 10 percent for over twelve 
hours. Occasional larger fluctuations seemed to 
be caused by local variations in the check side 
arm or check probe and to have no relation to 
the H concentration near the main side arm, as 
this probe did not show fluctuations at the same 
time. The best check on our temperature vari- 
ation curves however is their reversibility. In 
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runs extending even over several days, with 
furnace temperature raised from 20°C to 580°C, 
the probe temperature would drop to } its 
original value or less, and on returning the 
furnace to 20°C, the probe temperature would 
return to within 20 percent of its original value. 

From the known rate of loss of atoms to the 
walls at room temperature and the observed 
change in probe temperature with temperature 
of section L of the side arm, the variation in y 
may be found by a simple computation providing 
we assume we can piece together solutions of the 
form n= exp —By'x rather than complete 
solutions of the form 


n=A exp By'x+B exp —By'x. 


This assumption is valid however whenever 
we are sure that we can neglect the reflection of 
atoms beyond the length L, and can therefore 
attribute all changes of m to the change of decay 
along the heated tube. Actually the length L is 
long enough so that at room temperature the 
value of 1 falls to a small fraction of its original 
value at the far end of L, and as actually the 
recombination coefficient was in all cases ob- 
served to rise with temperature, this condition is 
even better fulfilled at higher temperatures, At 
temperatures below room temperature, y remains 
nearly constant, again introducing no appreciable 
error. We therefore derive 


t=k’ exp (—BiyiL) ‘exp [ —Bovol(x—L) ] 
where k’ is some constant 


to=k’ exp (—Boyo'x) 
whence 


to/ti=exp (Bi1y'L —Boyo'L), 
loge (to/ts) = Biyi'L —Boryo'L, 
BoyeL +log, (to /ty) . 
Ti= ( ) 
BiL 
to= probe T for L at To, 
t;=probe T for ZL at 7;. (18) 


TABLE I. Variation of 8 with temperature. 














rc B F di 8B 
— 180 7.26 180 2.91 
— 80 4.70 300 2.59 
15 3.73 450 2.28 
590 2.09 








The variation of 8 with temperature must also 
be computed. We remember 


B= (6C/4SD)}. (19) 


With p=75y, using Amdur’s" value for D, we 
obtain Table I for 6. 

The experimental runs were taken by varying 
the furnace temperature in steps of about 120°, 
allowing from one hour to five hours or more for 
the furnace to reach equilibrium. A minimum of 
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Fic. 5. Variation with pressure (in yu) in the ratio of 
temperatures of two probes at two different distances from 


a He discharge. Circled point corresponds to circled points 
of Fig. 4. 


one hour was found necessary for any approach 
to equilibrium, and several runs were made 
leaving the furnace at various temperatures 
overnight to be sure a close approach to equi- 
librium was attained. As the variation of probe 
temperature with furnace temperature was 
reversible, the values obtained do represent 
equilibrium values. Several different pressures 
near 75 microns were used, and since the results 
for all the different runs showed more random 
variations than systematic variations, all were 
averaged and plotted in Fig. 3. The low tem- 
perature portion of the experiment was con- 
ducted with the section L bent into a U tube and 
immersed in dry ice or liquid air. 


117, Amdur, J. Chem. Phys. 4, 339 (1936). 
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Fic. 6. 1, y(H) on Pyrex in H2+2 percent H:O at 75y 
pressure as a function of temperature. 2, y(OH) on Pyrex 
in H2O at 50y pressure as a function of temperature. 


The value of y may be calculated from Fig. 3, 
Eq. (18), and Table I. The results are plotted in 
Fig. 6, together with those for y (OH) on Pyrex 
(see Section 5). We may note that the observed 
variation cannot be simply explained by any 
mechanism involving some single energy of 
. activation or heat of reaction as the resulting 
equations will predict log y inversely propor- 
tional to temperature rather than almost directly 
proportional as observed in the high temperature 
region of Fig. 6. So simple an explanation could 
not be expected however as the temperature is 
sure to modify the properties of the surface, in 
particular the layer of adsorbed water. 


E. Variation of y with Water Vapor Content 


Three runs were made with respectively } 


percent, 2 percent, and 9 percent water vapor 
added to the hydrogen. It was found that the 
probe temperature increased with increasing 
water vapor content for all furnace temperatures, 
but especially at the higher temperatures. The 
results are plotted in Fig. 7. A pressure of 60 
microns was employed, and at least three hours 
were allowed for equilibrium at each temperature. 
The three curves were obtained in one con- 
tinuous run of over 200 hours; curves 2, 1, and 3 
consecutively. In such a long run, the probe is 
apt to change its sensitivity somewhat and there 
is reason to believe curve 2 should be more nearly 
intermediate between 1 and 3 than it is, as a 
check point taken at room temperature with 
H.+2 percent HO just before the run 3 gave 
a probe temperature of 175°. The general nature 
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of the variation with water vapor content is 
however unmistakable. 

Increasing water vapor content in the dis- 
charge may mean either better poisoning of the 
walls or increased OH concentration (see Part 5) 
or both. In Part 5 it will be shown that the 
uncertainty in OH decay is such that OH radicals 
might still reach the Pt probe. However, using 
the most probable value of OH decay, the OH 
concentration at the Pt probe may be neglected, 
and the observed variation may be attributed 
entirely to better poisoning of the glass walls for 
H recombination. 


4. RECOMBINATION OF H ON OTHER SUBSTANCES 


Substances investigated include quartz, com- 
bustion tubing, and K,SiO; chosen for their 
similarity to Pyrex; KCl because of its im- 
portance in certain investigations of OH recom- 
bination ; AlsO; as a typical dehydration catalyst ; 
ZnO-Cr.O3asa typical dehydrogenation catalyst; 
Pt; and also KOH, K2CO3;, metaphosphoric acid, 
and Na3PQx,. 

Two arrangements were used in the investiga- 
tion. The first is the same as was used in Section 1 
on Pyrex—an indirect method, good for substances 
of low catalytic power. In this method, the atoms 
diffuse down a side arm coated in part with sub- 
stances under investigation, and the atoms 
passing through this coated section are detected 
by a Pt coated probe (see Fig. 1). The coated 
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Fic. 7. Variation of Pt probe temperature with water 
vapor content of discharge. Arrow points to an independent 
check of first point of curve 2. 1=H2+9 percent H,0. 
2=H2+2 percent H.O. 3=H2+}3 percent H:0. 
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TABLE II. H recombination on various substances. We may note from the close check between the values of the last 
two rows that vy (Pt)200° is nearly equal to 1. 







































































Probe U 
temp. 
corrected 
Pt for 
Sample probe Probe U changing 
length L temp. temp. atomic ¥ v/v Pt 
probe Furnace (indirect (direct concentration (indirect (direct 
Substance dist. x Condition temp. method) method) and speed method) method) 
Quartz L=20 erratic 15 45 —- — 7.0-1074 
x= 30 behavior 300 22 a — 2.6-1073 _- 
500 15 -— — 3.4-1073 -— 
Combustion 20 15 70 oo — 4.6-10~4 _— 
Glass 30 300 50 -— — 1.8-10-3 — 
500 30 aoe = 2.8-1073 = 
KSiO; 6 moist 15 200 0 0 10-8 1.7-10-3 
30 dry 15 10° — —- 7-107 oa 
dry 200 4° 35 47 7-10 1.6-107" 
dry 500 i — 1-107! — 
KCI 6 15 350 0 0 2-10-75 
30 400 240 0 0 2-107 
KOH  }}}73=—” — moist 15 200 0 0 10>! 10-3 
30 dry 15 0 — = = -107 — 
dry 200 0 20 27 — 9-10 
dry 300 0 _ —- =>2-1071 — 
K2CO; 6 moist 15 200 0 0 10-4 10-3 
30 dry 15 0 14 14 =>1-107 5-10? 
dry 200 0 10 13 =>2-107 4-107? 
dry 400 0 — _ =>2-107! a 
Metaphosphoric acid 20 moist 15 350 0 0 2-10-75 os 
30 
Sodium 6 moist 15 200 7 0 - 10-1 =10-3 
Phosphate 30 moist 300 150 — — 5-10-4 — 
dry 15 0 20 20 =>1-10— 7-107? 
dry 200 0 15 20 — 7-107? 
dry 300 0 — — =2-107 — 
Al,O; 25 moist 15 20 20 20 9-102 7-107? 
30 moist 200 _ 10 13 — 4-107 
dry 15 0 100 100 - 3.3-107! 
dry 200 0 60 80 — 2.7-107 
dry 400 0 25 67 _— 2.2-107! 
ZnO-Cr.0; 2.5 moist 15 0 160 160 — 5.3-107 
30 or 200 0 120 160 — 5.3-107! 
dry 400 0 65 170 -= 5.6-107! 
Pt 2:5 moist 15 0 300 300 a 1 
30 or 200 0 250 330 — 1.1 
dry 400 0 120 320 - ia 











1 Numerous other runs established that within experimental error the catalytic effect of KCl] on H recombination cannot be distinguished from 
that of Pyrex. As the KCl coating did not completely shield the glass surface from the discharge, it is therefore possible that KCI is even less 
catalytic to H recombination than Pyrex. This is hard to establish definitely as measuring the recombination coefficient on glass represents about 
the limit of our experimental technique with the present method of investigation. 


sections were 11 mm inside diameter, glass blown 
to the system save for quartz, combustion tubing, 
Al,O3, ZnO-Cr.O3 and Pt. (These latter sub- 
stances were of 8 mm inside diameter tubing 
slipped inside the main side arm where they 


fitted snugly.) 


The indirect method is obviously not good for 


substances of high catalytic power such as 





K.CO;, KOH, etc., because few atoms would 
reach the probe, but we have used it as a very 
rough estimate of y. Where zero probe readings 
were observed, we have assumed the probe 7 to 
be equal to or less than one-half degree and we 
have assumed that the probe T at the discharge 
end of the tube would be at least 500°. These 


assumptions yield our lower limits of y listed in 
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Table II. y is calculated from the probe tem- 
peratures, the data for Pyrex, and the obvious 
equation” 


log. (t) = — (ByL)2—[By(X —L) Jpyrex +C’ 


where z is the substance under investigation. 

The second system of probes, utilizing a direct 
method of comparing different substances, is 
illustrated in Fig. 8. It is good for substances of 
high catalyzing power. 

The unknown substance is coated on probe 
“U.” Probe ‘‘T”’ is coated with Pt as an indicator 
of the minimum H concentration near ‘‘U.”’ (The 
actual concentration near ‘‘U”’ is of course in 
general much greater than that indicated by 
“T,” and the reading of ‘‘T’’ was never inter- 
preted quantitatively. In some runs however, 
a portion of the active substance coated on probe 
““U”’ would become dislodged or would be sub- 
limed onto the surface of the side arm. This 
coating then would decrease the total H concen- 
tration as measured by ‘‘7”’ without increasing 
the reading of ‘“‘U.’”’ The presence of the check 
probe “7” was useful in eliminating these runs.) 

Probe ‘‘U”’ is inserted in the middle of a 20 
cm long furnace whereby the temperature vari- 
ation in the response of ‘“‘U’’ can be measured. 
To correlate the readings of probe ‘‘U”’ at dif- 
ferent furnace temperatures, the change in heat 
loss of the probe with temperature of the sur- 
rounding furnace must be known. Calculations 
of heat losses, including radiation losses,— 
necessarily rough because of the comparatively 
long mean free path (of the same order of mag- 
nitude as the probe diameter) and the inac- 
curately known accommodation coefficient— 
show that this heat loss does not vary very 
greatly with temperature over the range from 
20°C to 450°C. That is, in Eq. (12), the term A/k 
may be taken as constant over this temperature 
range, and g will be proportional to y. In this 
equation however, the factor nd varies with tem- 
perature. We assume 1 is not different for dif- 
ferent coatings on the probe. This assumption is 
only a first approximation, and as n will actually 
be somewhat larger in the presence of less 
catalytic substances, a direct comparison of these 
with a Pt coating on “U”’ will make the unknown 


(20) 


122 Derived analogously to Eq. (14). 
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Fic. 8. Direct method of investigating y. Substance to be 
investigated is coated on ‘‘U”’ probe. 


substance appear more catalytic than it really is. 
N also varies with temperature however due to 
losses to the clean glass walls of the side arm. 
This correction may be computed from Eq. (20), 
substituting L=10, X=17 (here z is Pyrex at 
the furnace temperature), the observed tem- 
perature variation of y (Fig. 6), and the com- 
puted values of 6 (Table I). These changes in 
and @ are corrected for in Table II. 

The data in this section are not as accurate as 
in Section 1. In the first place, the salts at high 
temperatures (around 500°C) sublimed some- 
what to different parts of the side arm, and even 
a slight invisible layer of salt on the glass com- 
pletely changes its catalytic effect. In the second 
place, certain salts such as KCl were reduced in 
or very near the discharge. Also, some substances 
lost water of crystallization in violent miniature 
explosions which spread the coating. Most im- 
portant of all, however, was the fact that the 
catalytic effects of almost all these substances 
depended far more on their moisture content than 
on their temperature. Moist salts were non- 
catalytic; dry ones (save KCl) catalytic to H 
recombination. As a result, the data obtained 
with these salts with varying temperature are in 
general non-reversible, low catalytic effect being 
initially observed, increasing as the increasing 
furnace temperature dries out the salt, and 
remaining high when the furnace temperature is 
again lowered. 


5. EXPERIMENTS WITH WATER VAPOR 


When a discharge is passed in water vapor, the 
molecules are dissociated into H atoms and OH 





iC 











SURFACE RECOMBINATION OF RADICALS 121 


radicals. There is a negligible concentration of 
O atoms.5 Three possible two-body surface 
recombinations can occur. 


i— H+H+M—H.4+M/+4+4.3 ev,® 
2— H+OH+M—H,04+ M455 ev, 
3—OH+0H+ M—H.0.4+ M+3.1 ev. 


M is the surface catalyzing the recombination. 
All three reactions are strongly exothermic. No 
previous work has been done to distinguish 
between recombination processes 2 and 3. 

The state of knowledge with respect to these 
three recombinations prior to these experiments 
may be briefly summarized. All the data were 
known only at room temperature. At this tem- 
perature Pyrex glass was known to catalyze 1 
very poorly and 2 and 3 also poorly, although it 
was thought to be somewhat more catalytic for 
one or both of the latter reactions than for reac- 
tion 1. Numerous substances—such as Pt, 
ZnO-Cr203, KOH, etc., were known to catalyze 
1 strongly, but their effect for 2 or 3 was un- 
known. Al.O; was known to catalyze 1 strongly 
and was thought to catalyze 2 also. KCI was 
known to be quite ineffective in catalyzing 1, and 
quite effective in catalyzing either 2 or 3 or both. 
All these conclusions may be drawn from the 
work of Taylor and Lavin,’ although they do not 
try to distinguish the two possible types of OH 
recombination. 

The present investigation has added to the 
existing data and narrowed down the possible 
interpretations, but has not succeeded in 
uniquely distinguishing between recombination 
processes 2 and 3. We have succeeded in obtain- 
ing a rough estimate of the recombination coef- 
ficient y for process 2 or 3 on glass, and have 
found its change with temperature. We have 
also investigated the effect of temperature on 
other substances, especially KCl. 

The experimental arrangement is indicated in 
Fig. 1. Water vapor is pumped through the dis- 
charge tube from the water reservoir. The OH 
concentration is detected by the heating of 
probes coated with KCI. (Note that Pt probes 
could not be used as Pt recombines atomic H. 
KCI does not recombine H and does recombine 


8 Energies of recombination calculated from data of 
Int. Crit. Tab. V, p. 176, and A. v. Engle and M. Steenbeck, 
Electrische Gasentladungen I, p. 43. 


OH—see below.) The response of these KCl 
surfaces was quite erratic, and they almost always 
required activation for 15 minutes at 250°C to 
manifest their catalytic power. The probe 
response was fairly reproducible for an average 
of about 5 runs, after which it would usually 
become very erratic. The glass walls near the 
probe were observed to become coated with KCl 
in patches by fragments exploded from the 
probe in the course of repeated activation. This 
coated section robbed the probe of OH and 
caused part at least of the erratic nature of the 
probe response. Because of this effect, probe 
temperatures below 10° were considered unre- 
liable. 

To avoid possible reduction of the KCI by the 
discharge, the probe was always kept at least 10 
cm distant from the discharge. 


A. Pyrex 


1. Investigation of y (OH) at 
Room Temperature 


Figure 9 is a logarithmic plot of KCI probe 
temperature against distance from the discharge. 
The scattering in observed points is very large. A 
greater range of distance could not be investi- 
gated as for shorter distances there was danger of 
reduction of KCI by the discharge and for larger 
distances the probe temperature was so small as 
to be unreliable. The scattering here, in contrast 
to that in Fig. 4, is believed due more to variation 
in probe response than in surface conditions of 
the glass. One sequence of points (the circled 
values) was taken with the same KCI coating on 
the probe, readings at small and large values of 
x alternating. 

Whether OH disappears by process 2 or 3 is 
immaterial in the calculation of y as y is defined 
as the fraction of times the OH radical strikes 
the wall and is recombined there by any process. 
We have therefore the same equation as for 
7y(H). [Eq. (14) ], but in this case Bp =4.63 (see 
Section 2). From Fig. 9, we get therefore 


y=8X10-5. (21) 


The small range in probe temperatures (a 
factor of 10) and distances (a factor of 5) avail- 
able (compare factors of several hundred and 20 
respectively in experiments on H recombination) 
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Fic. 9. KCl probe temperature as a function of distance 
from a water vapor discharge. Pressure=50yu. Circled 
points are for one sequence of runs with the same KCl 
surface on probe. 


together with the relative unreliability of the 
KCI probes rendered attempts to distinguish 
between surface recombination and losses of OH 
by three-body collisions useless. However as the 
pressures are even lower than in the work on H 
recombination, surface effects probably pre- 
ponderate. The nature of the temperature vari- 
ation in probe response (Fig. 10) evidently 
precludes the possibility of recombination at the 
probe overshadowing losses of radicals to the 
Pyrex walls. 


2. y(OH) at Other Temperatures 


Investigations were made in the temperature 
range from 15°C to 590°C. The experimental 
set-up is the same as that in Fig. 1 except that 
two probes ‘“‘U,”’ and “U,”’ in two independent 
side tubes were used instead of one. The double 
probe system was used to eliminate random 
fluctuations. An average of several runs is plotted 
in Fig. 10. Here the temperature of the 7 cm 
furnace was varied by jumps of 150° every two 
hours. Reversibility of probe response with 
furnace temperature showed that probe readings 
represented equilibrium values. The short reliable 
life of the KCI probes made it impossible to vary 
the furnace temperature more slowly. 

Computation of 6 for OH diffusion is not ac- 
curate as no data are available on the coefficient 
of diffusion of OH. We may get a rough estimate 
however by assuming that the coefficient of 
diffusion of OH nearly equals the self-diffusion 
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coefficient of H,O. This we compute from the 
viscosity data by means of the equation 


D=Cn/p (22) 


where we take C= 1.36 as a compromise between 
the hard sphere and inverse fifth power repulsion 
values of C (1.200 and 1.543, respectively). The 
resulting values of 8, calculated from equations 
analogous to those of Part 3, are givenin Table III. 

The values of y computed from this table and 
the data of Fig. 10 are plotted in Fig. 6. 


3. y(H) ‘on Pyrex in HO 


In Part 3 we have measured y(H) on Pyrex 
in moist Hz. Here we may do the same for H 
atoms produced in a water vapor discharge. Dif- 
ferent values of coefficient of diffusion, etc., make 
this determination of y an interesting check with 
the previous ones. Furthermore, if the recom- 
bination H+OH—H:0O is considerably stronger 
than either H+H—-H,. or OH+OH—H.O,, a 
large value of y should result as the initial dis- 
sociation of H,O yields equal amounts of H and 
OH. 

OH disappears more rapidly with distance 
from the discharge along the side arm than does 
H. This is known from the data of Taylor and 
Lavin,’ and also from the large value of y(OH) 
which we have found. It is not known whether 
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Fic. 10. KCl probe temperature as a function of furnace 
temperature in a water vapor discharge. Pressure = 40x. 














ce 








SURFACE 


OH recombines on Pt or not, but presumably it 
does. Consequently the Pt probe is affected by 
both H and OH. However, for distances suffi- 
ciently far from the discharge (around 50 cm), 
the larger recombination of OH on Pyrex assures 
small concentrations of OH compared to H, and 
thus the decay of H alone may be measured. 

Pt probe readings in H.O are therefore plotted 
logarithmically in Fig. 11 for a pressure of 50x. 
The scattering of points is not great. The devi- 
ation from a straight line shown by the initial 
group of points is presumably from the effect of 
OH on the probe. The final two points near 300 
cm are zero probe readings whose logarithms are 
plotted arbitrarily as zero (equivalent to a tem- 
perature of 1°), and their deviation from linearity 
is not significant. 

The following equation" holds for the diffusion 
of molecular He into HO. 


D=0.7516(T/273)p/160 (pinmmofHg). (23) 


As H has nearly the same diameter as H, © 
and as for large differences of masses of the inter- 
diffusing molecules, D is proportional to the 
velocity of the lighter component, we can con- 
clude that for the diffusion of atomic H into H2O, 


D=2}-0.7516(T'/273)p/760. (24) 


With this value of D, for a pressure of 50, at 
room temperature, we get 8=3.72, whence, using 
Fig. 11 and Eq. (14), we get 


y(H) =1.9X10-. (25) 


The agreement with values of y for moist 
hydrogen is surprisingly good and indicates that 
the recombination H-+-OH—-H:0 on Pyrex 
is weaker than H+H—H, or OH+OH—H.0.. 

Experiments similar to those in Part 3 showed 
three-body collisions to be negligible compared 
to the surface recombination of H. 


TABLE III. Values of 8 for OH diffusion. 











T°C 8 T°C r) 

20 4.63 450 2.08 
180 3.16 590 1.89 
300 2.53 








4 TInt. Crit. Tab. V, p. 62. 
% Calculations from Amdur’s viscosity data (reference 
11) for atomic H give this result. 
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B. Other Substances 
1. KCl 


Numerous experiments with alternating He» 
and H,0 discharges establish definitely that KCI 
coated probes are heated appreciably—up to 
around 100°—at 10 cm from a H.O discharge 
and are not heated above 1°-2° at the same 
distance from a discharge in moist He. We may 
therefore conclude that KCl catalyzes OH 
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Fic. 11. Pt probe temperature as a function of distance 


from a water vapor discharge. Pressure=50. O probe 
temperatures plotted as 1°. 


recombination and does not catalyze H recom- 
bination. This is important because it establishes 
a test for free OH which is not affected by H. 
Our conclusion agrees with the results of Taylor 
and Lavin’ and Frost and Oldenberg.'® 

Similar experiments with a furnace surround- 
ing the probe have been carried out. At a furnace 
temperature of 180°C, the KCI response de- 
creased by a factor of two. Possibly this was 
due to some KCI sputtered onto the walls. Inser- 
tion of a check probe beyond the furnace proved 
useless as the OH concentration was too small 
to be accurately detected. 

To estimate the magnitude of y(OH) on KCl, 
we observe that a KCl coated probe 20 cm from 
a H.O discharge heats up to 65° whereas a Pt 


16 A, A. Frost and O. Oldenberg, J. Chem. Phys. 4, 642 
(1936). 
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TABLE IV. Comparative probe temperatures with 
He and H.0O discharges. 








AlsOs AleO3 
Gas (unactivated) (activated) 





ZnO -Cr203 Pt 
He 290 315 290 580 
H,O 80 180 220 415 








probe 20 cm from a He: discharge heats up to 
500°. As y(OH) on Pyrex is greater than y(H) 
on Pyrex, the OH concentration near a HO 
discharge surely does not exceed the H concen- 
tration near a He discharge. Furthermore the 
mean of the heats of recombination by the 
processes H+OH—H.,0 and OH+OH—H.0O, 
is nearly the heat of recombination of atomic H. 
We may therefore conclude that y(OH) on KCl 
is surely greater than 0.1 times y(H) on Pt. 

The temperature variation of OH recombina- 
tion on KCl could not be followed closely, but 
since a KCI probe located 10 cm beyond a 6 cm 
section of the side arm coated with KCI showed 
no response for any temperature of the 6 cm 
section up to and including 400°C, evidently 
KCI remains catalytic to OH recombination at 
high temperatures. 


2. KOH, K;COs3, K.SiOs;, Nas3PQO, 


These substances all have the property of 
recombining atomic H strongly when dry and 
poorly when wet. Measurements of OH recom- 
bination by the direct method therefore are useless 
as the H recombination masks the OH response, 
and any systematic variations between response 
in He and in H.O are as easily ascribed to varying 
catalytic effect for H (because of the strong 
dependence of H catalysis on surface dryness) 
as to some catalytic recombination of OH. 

Measurements by the indirect method with KCI 
coated probes are useless if the KCI catalyzed 
recombination is by the process H+OH—H,0O, 
as no H passes the coated section of the side arm. 
Similarly, measurements with a Pt probe are 
valid only if the reaction OH+OH—-H,0, 
proceeds strongly on Pt. The nature of OH 
recombination on KCl and Pt is not known 
however. Furthermore, a positive response with 
a Pt probe can also be interpreted as a change in 
7y(H) on the coatings on the Pyrex walls with 
increased H,O content of the discharge. 


SMITH 


We have nonetheless made these measure- 
ments, as a positive response of the KCI probe 
in any case where no H passes the coated section 
would establish the reaction OH +OH—H,20: on 
KCl, and a null response with the Pt probe 
shows y(H) is still large on the wall coatings 
investigated. 

The substance under investigation was coated 
inside the 7 cm length of side arm in the furnace. 
Probes ‘‘U,’”’ and “U2” were located 5 cm beyond 
the furnace, and were coated with KCI. Probe 
“T’’ was located 5 cm beyond U; and U2 and 
was coated with Pt. Alternate runs were taken 
with a clean glass side arm to check the response 
of the KCI. 

Positive readings on U; or U2 were found only 
with a NasPQ, coated side arm at temperatures 
below 400°C—i.e., in an unactivated condition 
where H_ recombination was known to be 
uncatalyzed. This therefore established only the 
fact that unactivated Na3PQ, catalyzes neither 
H recombination nor OH recombination. 

Null readings were found with “7” for 
activated coatings, except that readings of a few 
degrees were observed with a K,SiO; coated side 
arm even at high furnace temperatures. This 
behavior is the same as with moist hydrogen, 
showing that for the substances investigated it 
is a function of surface wetness, not H2O concen- 
tration in the gas. 


3. Pt, AlsO3 and ZnO-Cr.O3 


Indirect measurements with these substances 
(3 cm samples) established zero KCl probe 
readings which proved little (see previous sec- 
tion). The Pt probe response was the same as 
with Hz (part 4) except that with an Al,O; 
sample, Pt probe readings of 1°-2° were observed 
at room temperature even after activation of the 
Al,O; at high temperatures (at which tempera- 
tures the Pt probe response was zero). Probable 
interpretation of this result is that there is 
greater poisoning of Al.O; for H recombination 
in HO than in moist He at the pressures em- 
ployed. (See below for confirmation of this 
interpretation). 

Direct measurements at .17 cm from the dis- 
charge established the data in Table IV for alter- 
nate H, and H.O discharges. (The data are aver- 
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age values for several runs. He pressure 754; H2O 
pressure 50x). 

To obtain high readings with an H,O dis- 
charge, it was necessary to heat the Al2O3 probe 
in a furnace to 320°C for 30 minutes. As even 
after this activation the Al,O; probe temperature 
slowly decreased over the course of 30 minutes 
to 125°, it seems probable that the lower ratio 
of readings with H.2O to readings with Hy» for 
Al,Os compared with ZnO-Cr.O3; and Pt is 
caused by increased poisoning of the AIO; 
surface by water vapor. 

Our observation that Al,O; does not heat up 
more with a H2O discharge than with a H, dis- 
charge disagrees with that of Taylor and Lavin. 
However if the ratio of their Al,O; probe readings 


to their silver probe readings is examined, it is 
seen not to change appreciably with changing 
H;0 concentration in the discharge. Their data, 
together with ours, therefore show no significant 
difference in relative heating of AloO;, ZnO-Cr2QOs, 
Pt, and Ag with H.O and He discharges—and 
therefore no measurable difference in the relative 
catalytic effect of the dehydration catalyst AlsO; 
and the dehydrogenation catalyst ZnO-Cr2O; for 
H and OH recombination. 

The author wishes to express his appreciation 
to Professor O. Oldenberg for suggesting this 
problem and for his continued interest and help 
in its development. He is also indebted to Pro- 
fessor W. H. Furry for clarification of many 
aspects of the diffusion problems considered. 
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A molecular model in terms of which the elastic viscous properties of rubber-like substances 
can be interpreted is presented. Experiments on stress relaxation at constant extension, creep 
under constant load, extrusion, vibration, and breaking are discussed in terms of a mathe- 


matical formulation of this molecular model. 


VERY large class of substances composed 

of large molecules, such as nylon and other 
fibers, natural and synthetic rubbers, muscle 
tissue, and so on, exhibit elastic-viscous proper- 
ties that are qualitatively similar, which fact 
reflects an over-all similarity in molecular 
structure. 

These substances are composed of long chain 
molecules linked together at some points by 
primary cross bonds and at other points by 
bonds of a secondary nature such as strong 
Van der Waals’ bonds (as in regions of local 
crystallinity) dipole-dipole bonds, etc. The cross 
linking bonds give these polymers an effective 
network structure. Between the cross bonds 
there is a rather long length of coiling, twisting 
molecular chain, the segments of which are 
rotating more or less freely. This amorphous 
liquid-like structure comprises the bulk of the 





polymer, when the latter is in the temperature 
range of the rubbery state of aggregation. 

Because of statistical considerations, the 
molecular chains spend most of their time in very 
coiled configurations when no external stresses 
are acting. When the polymeric network is 
rapidly stretched, the molecular chains are 
stretched and thereby have fewer conformations 
that are consonant with this extended state. 
The probability of a unit cell of the network 
having dimensions x, y, 2 is'? 


Prey) =A exp [—B(x?+y?+2") J. (1) 
The entropy is related to the probability by 


1E. Guth and H. M. James, Ind. Chem. 33, 624 (1941). 
See also earlier work of Meyer, Mark, Kuhn, and Guth 
referred to in this paper. 

2 Tobolsky, Powell, and Eyring, Chapter on “Elastic 
Viscous Properties of Matter” in Frontiers of Chemistry, 
Interscience Publishers (in press). 
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Boltzmann’s relation S=kinP, and the tension X 
to the entropy by the relation L= —7(0S/dx) r. 
It is assumed that the volume of the unit cell of 
the network xyz=v does not change during 
stretching because of the interpenetrating char- 
acter of the network structure. On this basis it 
had been shown by Guth and James! that the 
tension X is given by 


Pn ccicansil L =I" (2) 


Lo '!Lo (tL 
Choosing suitable values of Z, the number of 
unit cells in the network, one can explain the 
values of the modulus of rubber which is about 
10’? dynes/cm? instead of 10" dynes/cm? char- 
acteristic of short chain crystalline substances. 

Many of the interesting phenomena regarding 
the elastic behavior of rubber are described by 
such terms as creep, creep recovery, relaxation, 
hysteresis, etc. We shall treat these phenomena 
in terms of molecular rate and flow processes. 

In any given cross section of the polymer 
perpendicular to the stress there are a number of 
deformation mechanisms in which relaxation 
processes can occur independently, and over 
which the stress distributes itself. In all these 
unit processes of deformation, the effect of 
stress on the rate of strain is twofold, there 
being first an elastic component for which the 
rate of strain is in phase with the rate of stress, 
and a ‘contribution to the rate of strain of 
relaxation processes such as slipping, tearing, 
etc.? It is convenient to classify the various unit 
processes of deformation into categories which 
group together similar processes. Let us suppose 
that there are in a square cm of cross section N; 
such parallel processes of type i, so that if the 
total stress f; distributes itself on this type of 
unit process of deformation, the average force 
acting on each is f;/N;. Suppose, too, that per 
unit length, there are m; such processes of 
deformation in series, then the total rate of strain 
is given by 








ds 1 df; kT 
me hy 
dt G; dt h 
iNi 
Xexp (—AFF/kT)2 sinh , & 
Pe / 2NikT 
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where \; is the average distance projected in the 
direction of stress between equilibrium positions 
in the relaxation process and AF* the free energy 
of activation for the process.* G; is the appro- 
priate partial modulus to be associated with the 
deformation process, and the first term on the 
right represents the elastic contribution. The 
contribution to rate of strain due to relaxation 
(slipping, flow, etc.) is given in the second term, 
and this term is derived in similar fashion to the 
treatment of viscous flow. The force f;/N; acts 
on each unit of type 7 and tends to favor motion 
from one equilibrium position to another in the 
direction of stress. The sum of all the partial 
stresses on the various deformation mechanisms 
is equal to the total stress and the stress is 
distributed in such a way that the rate of 
strain is the same for each. 

Equation (3) can be written in a more con- 
densed form: 





ds 1 df; 
—=— —+4A,sinh B,f,. (4) 
dt G; dt 
In case B;fsK1 then 
ds 1 df; 1 df; 1 
a —$hRi——— ef, ) 
dt G; dt G; dt; 


Equation (5) is of the form that was postulated 
by Maxwell to explain elasto-viscous behavior. 

The deformation mechanisms found in poly- 
mers in the rubbery state of aggregation are of 
the following types. There are first the structural 
elements of the network, which support most of 
the stress. Some of these elements are composed 
of primary cross-linking bonds, and so behave 
purely elastically for small loads. Other units of 
the network structure are partly joined together 
by secondary cross bonds, which can slip on the 
application of a stress. Some of the total stress 
is also distributed on the mobile segments of the 
long molecules which form an amorphous liquid- 
like structure pervading the network. The 
equations of motion for the polymeric material 


3 The flow term in (3) will be seen to be the same as (8) 
in a previous paper [H. Eyring, J. Chem. Phys. 4, 283 
(1936) ] if one writes m=),~!, Ni; =(AzA3), and A;=A. 
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are: 
ds 1 dfi : 
(a) —=— —+A,sinh B;f, 
dt G, dt 
primary network units, 
$a hang (6) 
»)) = — sin 6 
dt Godt = 
secondary network units, 
ds 1df; 1 
(c) —=— —+—fs segment motion. 
dt G3; dt 73 


The constants A and B are to be interpreted 
in terms of Eq. (3). A; is small, and in some 
cases the breaking of primary cross bonds is to 
be neglected. G,; and G: are moduli of the 
“entropy elasticity” type discussed previously, 
and of the order of 10°—108 dynes/cm?. G, and 
Ge are related to m;, Ni, m2, and N..* G; is the 
same kind of modulus as is associated with 
short chain crystalline hydrocarbons, namely, 
of the order of 10" dynes/cm?. The viscosity 
associated with the motion of segments is much 
smaller than that associated with slipping at 
secondary network cross bonds, and since the 
stress distributed on this deformation mechanism 
is small, the equation is of the Maxwell type. 

A scheme mathematically similar to the above 
had been introduced by Kuhn‘ and extended by 
Simha.’ However, these authors use Maxwell 
type equations throughout. A non-mathematical 
discussion of the properties of hard rubber has 
recently appeared which uses a qualitative 
picture and terminology similar to that used 
above.® A somewhat similar qualitative picture 
has been presented by other workers for various 
polymeric materials.” § 


t m2 is the reciprocal of 1, the mean distance between 
secondary network bonds projected in the direction of 
stress. Ne is the number of secondary bonds per unit area 
in a slab whose height is /. The number of secondary bonds 
v2 per unit volume is equal to m2N2. Similar considerations 
hold for primary bonds. 

*W. Kuhn, Zeits. f. physik. Chemie B42, 1 (1939; also 
K. Bennewitz and H. Rotger, Physik. Zeits. 40, 416 (1939). 

5R. Simha, J. App. Phys, 12, 680 (1941). This paper 
gives a large list of references to early work in this field. See 
also J. D. Ferry, J. Am. Chem. Soc. 64, 1330 (1942). 

6]. R. Scott, Trans. Faraday Soc. 38, 284 (1942). 

7 T. Alfrey and H. Mark, J. Phys. Chem. 46, 112 (1942); 
H. Mark, Ind. Eng. Chem. 34, 449 (1942); Mark, Cold 
Tout Harbor Symposia on Quantitative Biology 9, 204 

8 H. Leaderman, Lecture presented at A.C.S. meeting in 
Buffalo, New York, September, 1942. 


Equations (6), when applied to experimental 
results, present the possibility of analyzing data 
in terms of molecular properties. For example, 
B, depends on Ne, the number of secondary 
bonds per unit area. This will be seen to be a 
function of temperature and of initial deforma- 
tion. In the subsequent paragraphs we shall 
discuss some typical experiments which can be 
interpreted in the light of Eqs. (6). It is to be 
realized that very strong secondary bonds can 
play the role of primary bonds in certain 
experiments. 


STRESS RELAXATION OF CONSTANT 
ELONGATION 


If a substance is rapidly stretched to a given 
elongation and held fixed, the stress necessary 
to maintain that elongation decreases with time. 
For slightly vulcanized rubber Meyer and Ferri® 
concluded that the tension decays to a certain 
value below which the decrease is very slow. 
The same thing was observed by Misch and 
Picken’® for polyvinyl acetate. As a result of 
their studies of stress relaxation and decay of 
birefringence of polyvinyl acetate held at con- 
stant elongation these authors conclude that 
the polymer is composed of a netted system of 
chains through which interpenetrates a system 
of relatively free chains. 

For this type of experiment ds/dt=0, the 
stress due to the segment motion [Eq. (6c) ] 
decays so rapidly that it is not observable. The 
stress on the primary network units decays very 
slowly, and if a limiting stress is really reached it 
must be due to the primary network. We shall 
assume, therefore, that the primary bonds are 
not broken during the course of the experiment. 

The observed decay of stress is, therefore, due 
to the relaxation of secondary network bonds. 
If Maxwell’s equation held, the decay of fz would 
be given by f2=f2%e~“/", where 7 is the relaxation 
time (r=72/G2). The decay of stress does not 
follow such a law. Of course one can by using 
the Laplace transform find a suitable distribution 
of relaxation times, but it will be shown that 
integration of Eq. (6b), which has a natural 


9K. H. Meyer and C. Ferri, Helv. Chem. Acta 18, 574 
(1935). 

1 L. Misch and L. Picken, Zeits. f. physik. Chemie B36, 
400 (1937). 
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basis in terms of molecular behavior, can also 
explain the observed data. The result of such an 
integration is 


tanh (Bof2/2) 
=tanh (Bof.°/2) exp (—A2BeoGot), (7) 


where f.° is the initial stress. If Bof2/2 is larger 
than unity and f,°>f>1/2B2, one gets the 
approximate result that 


fo= —(1/B2) loge (¢2/2) —(1/Bs) loget, (8) 


where A2B2G2 is put equal to c. The total 
stress is, of course, fi+f2 and the decay of f; is 
very slow. 

It had been observed by Phillips'! that the 
decay of tension at constant extension for rubber 
can be represented by a formula of the type 


f=a—5b logt (9) 


in agreement with (6). 

In Fig. 1 is shown a plot of f2 against logo ct 
according to (5) for various values of the 
parameter B2f.°. Some of the subsequent figures 
show data for some of the stress relaxation 
measurements found in the literature plotted 
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Fic. 1. Relaxation of stress. 





uP, Phillips, Proc. Phys. Soc. London 19, 491 (1905). 








TOBOLSKY AND H. EYRING 


against logarithmic time. See, for example, Figs. 2 
and 3,!% 8 

By measuring the slope of the stress plotted 
against log time curves one can find values for Ne, 
the number of secondary bonds per unit area: 
this number is found to increase with decreasing 
temperature and also appears to be roughly 
proportional to the initial deformation or initial 
stress. 

The equation for stress relaxation 


tanh (Bf/2)=tanh (Bf°/2)e~* (10) 


appears to hold more generally than for poly- 
mers. In Fig. 3 is plotted data from Trumpler™ 
on relaxation of steel. The solid curves are 
calculated curves based on Eq. (10). The 
parameters B and ¢ for various values of the 
initial stress were very simply related. Bfo is 
found to be a constant (=6.32), and logc isa 
linear function of fo, which means that AFF is a 
linearly decreasing function of fo. Thus it appears 
that the number of relaxing mechanisms per 
unit area increases with increasing initial stress 
(or deformation) and the free energy of activa- 
tion for flow decreases linearly.* 


Creep under Constant Stress 


Several experiments in the literature deal with 
the creep of fibers or rubbers under constant 
load.* See Fig. 4. If the amount of creep is not 
too large, so that the change of cross-sectional 
area is not very great, we can suppose the 
experiments carried: out under conditions of 
constant stress. For small loads many of these 
substances show a delayed elasticity. There is 
an instantaneous deflection, then a recoverable 
creep till a limiting deflection is reached. Of 
course even for small loads, the primary network 
bonds will eventually break and flow. For nylon 
and other fibers for small loads and ordinary 
experimental times, a definite limiting deflection 
is reached. (See Fig. 5.) For natural and synthetic 


2L. R. G. Treloar, Trans. Faraday Soc. 36, 538 (1940). 

13. W. E. Trumpler, Jr., J. App. Phys. 12, 248 (1941). 

* This fact indicates that when the initial deformation is 
increased the units of flow become smaller (hence more 
numerous per unit area), and the free energy of activation 
decreases with the decrease in size in a manner to be 
expected. 

144 P, Kobeko, E. Kuvshinskij, and G. Gurevitch, Tech. 
Phys. U.S.S.R. 4, 622 (1937). 
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Fic. 2. Relaxation of stress at 25° for raw rubber (after Treloar). 
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remain rigid, and carry through our calculations 
on that basis. After a fairly long time the 
primary network units will also start breaking. 

For this case, too, the equation describing 
segment motion can be neglected. For a constant 
stress the immediate extension is seen to be 
f/(Git+G2). The limiting extension is f/G,, if 
breaking of the primary network is disregarded. 
We have to integrate the equations 


ds 1 dfi 
(a) ee ee ee 
dt G, dt 
6! 
ds 1 df. ( 
(b) —=— —+A? sinh Bofs, 
dt Ge dt 


under the condition that fi:+f2z=f, where f is 
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constant so that 


df;/dt= —df2/dt. 


Hence 
dfoy1 1 
a ee —-) =—As sinh Bafs, 
dt\G, Ge 
Let 
1 1 1 





tem, 
& G& G 


Integrating we get 











Bofs Bofe° 
tanh ——=tanh exp (—G,A2Bot), (11) 
where 
. Go f 
” GetGe 
From (6!a) 
S=(f—fe)/Gi. (12) 


The similarity of Eq. (11) to Eq. (7) is to be 
noted. 





For , 
fo>f>— 
~ ie 
1 1 
fro= sata log, A2B.G,— — log. ‘. (13) 
2 By 
so that 
f 1 
S=—+ log A2BoG,+—— log t. (14) 
Gi 2071 B 1 


It has been observed that the creep under 
constant load can be expressed in the form 
S=a'+b! log t. Once again we see that By is 
inversely proportional to initial stress, because 
b' above is found to be proportional to the stress. 

Figure 5 shows the creep under constant load 
on the basis of this theory for several values of 
the parameter Bf,°. It is to be noted that these 
values are closely related to those in Fig. 1 by 
reflection about the x axis. Since G,<G»2 the 
time scale for these types of experiments will 
be somewhat different. 

Data on natural ana synthetic rubbers creep- 
ing under fairly heavy load often appear to be 
best fitted by a sequence of straight lines when 
plotted against log time. This indicates that the 
bonds which have been assumed not to break 
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are also breaking after sufficiently long times. 
As before stated, there is probably a wide 
distribution of bond strengths, and the use of 
only two types of bonds is an idealization. 

The phenomena of creep recovery are best 
explained by the Boltzmann memory principle,'® 
which is a consequence of these equations of 
motion. 


Extrusion Experiments 


The non-Newtonian character of the flow in 
rubber-like materials is best shown by their 
behavior in being forced through a tube.? If the 
flow were Newtonian, Poiseuille’s equation would 
apply and the rate of efflux through the tube 
would be proportional to the pressure head. 
This is definitely not the case for plastic and 
rubber-like polymers. 

Suppose that we had a tube of length /, 
radius R, and were forcing material through it 
under a pressure head p. At any point whose 
distance from the axis of the tube is 7, the shear 
force f is given by 


f= pr/2l, (15) 


the shearing velocity du/dr we shall assume to 
be given by 


du/dr=A sinh Bf=A sinh (Bpr/2l). (16) 


By integration we can find the velocity in the 
axial direction as a function of r. We make use 
of the condition that the velocity is zero when 
r=R. This gives 


2Al BpR Bpr 
u == | cost —-—cosh — =| 
Bp 2/1 


(17) 


The volume of matter delivered per second is 
dV R 
—= f u2nrdr (18) 
dt 0 

and substituting for « from Eq. (17) 

dV “| BpR 





—= cosh —— 
dt Bp 21 


4l BpR 16/? 
_ sinh 
BpR 21 (BpR)? 


1% H. Leaderman, Textile Research 11, 171 (1941). 








BpR 
sinh? —| (19) 
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Fic. 5. Creep under constant stress. 


For sufficiently large pressures 


si son (BpR/2l) (20) 
—=—— exp (DPK/ 21). 0 
dt Bp 





Once again it must be remarked that whereas 
some of the data on extrusion indicate that 
Eq. (19) is essentially correct (see Fig. 6), other 
data for materials with a high percent of primary 
cross bonds or with large amounts of filler 
indicate that more than one type of flow 
mechanism has to be accounted for. The data 
in reference 16 appear to be fitted by N10" 
and AF*¥=15 kcals. 


Vibration Experiments 


If a rubber sample is subjected to a periodic 
stress, heat is generated because of the internal 
viscosities. Whereas the ‘“‘viscosities’’ that one 
calculates for experiments of the type described 
in the above paragraphs is about 10° poises, the 
viscosity calculated in vibration experiments is 
from 10*-10° poises. Most experiments on the 
effect of an alternating stress on rubber are 
carried out by a resonance method.'*'!8 A mass 
is attached to the rubber sample, and either the 
applied frequency or the mass is changed until 
the amplitude is a maximum, at which point the 
natural frequency of the system is equal to the 


16 J. H. Dillon and N. Johnston, J. App. Phys. 4, 225 
(1933). 

17S, D. Gehman, J. App. Phys. 13, 402 (1942). 

18 R. B. Stambaugh, Ind. Eng. Chem. 34, 1362 (1942). 
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externally applied frequency. From the height 
and breadth of the resonance curves one can 
calculate the ‘dynamic’ viscosity and ‘‘dy- 
namic’’ modulus of the rubber sample used. 

For these experiments the amplitude of 
deformation is quite small. The frequency is 
usually of the order of 60 cycle/sec. It is reason- 
able to suppose that large scale breaking of 
primary or secondary bonds is not involved, but 
that the viscous dissipation in these experiments 
is in large part due to a motion of chain segments. 

Another method of studying the relaxation 
mechanisms in polymers is by dielectric constant 
measurements of absorption and dispersion. 
Many: rubber-like polymers show a_ broad 
distribution of relaxation times of the order of 
10-® second. It seems unreasonable to attribute 
these to motions of entire molecules, and so 
these rapid relaxations must be due to motions 
of chain segments.!9 

Let us for simplicity assume that Kuhn’s 
formulation is adequate for these low stresses. 
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19W. Kauzmann and H. Eyring, J. Am. Chem. Soc. 62, 
3113 (1940). 
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Then in terms of the relaxation mechanisms 
the measured modulus and viscosity would be: 


iw? ;? 
G dynamically measured = }>; ————,,__ (21) 
1+w?r,? 
" Ni 
n dynamically measured = }); ———, (22) 
W°T 5" 


where w is the angular frequency and the sub- 
script 7 refers to the ith relaxation mechanism. 
We see from the above that 7 is a decreasing 
function of w and G is an increasing function of w, 
so that the fact that the “‘static’’ modulus is 
lower than the dynamic modulus measured by 
these experiments is not surprising. Also at very 
high frequencies we might expect that the 
moduli of rubber approach the values encoun- 
tered in crystalline hydrocarbons. 

If segment motion causes the heat dissipation 
in rubber subjected to rapidly varying periodic 
forces, we might then expect some kind of 
correlation between the effect of side chains on 
the internal viscosity of polymers and of the 
viscosity corresponding short chain compounds. 
Figure 7 shows the effect of adding various 
groups to a hydrocarbon skeleton six atoms long. 

Side groups alone do not determine the 
internal viscosity of polymeric materials. Inter- 
molecular forces play an important role, and if 
the polymeric chains, for stereochemical reasons, 
pack very well the internal friction is high. 

The exact dependence of G and 7 on frequency 
has not as yet been determined over a sufficiently 
wide frequency range to determine an appro- 
priate distribution of relaxation times for the 
segment motion. The indications are that the 
distribution is quite broad. 

If one makes the approximation that there are 
three relaxing mechanisms of the Kuhn type 
corresponding to the motion of primary network 
units, secondary network units, and segments, 
with relaxation times 71, 72, and 73, respectively, 
then Eqs. (21) and (22) can be simplified. Since 
71 is very large, 72 is of the order of seconds and 
73 is of the order of 10~* second, and remembering 
the values of the G’s and n’s discussed previously, 
it can be seen from Eqs. (21) and (22) that for 
frequencies of the order of 60 cycles per second, 
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G and 7 are approximately given by 
G=G,+G2=G dynamically measured, (23) 


Ne 


n ~-+n3 = dynamically measured. (24) 


1 +e?rs? 


In other words in a suitable frequency (perhaps 
from 20 to 200 cycles per second) the modulus 
should be independent of frequency and the 
viscosity plotted against 1/w? should be approxi- 
mately linear. Such appears to be the case for 
the data of Stambaugh.'* (See Fig. 8.) From 
the graphs discussed above one can find Gi;+Gz:, 
which is the “dynamic modulus.” 73 can be 
obtained as the intercept of the 7 vs. 1/w? plots, 
and G2/r2 is the slope of these plots. 

From Stambaugh’s data one obtains Table I. 
Since G; is approximately the “static modulus”’ 
as measured in non-alternating experiments, 
which is cf the same order of magnitude as the 
‘dynamic modulus” (often one-third to one-half 
as large), G2 can presumably be found. 72 will 
therefore turn out to be of the order of seconds. 

It is to be noted that both segment viscosity 
and viscosity due to relaxation of secondary 
bonds contribute to the dynamic viscosity at 60 
cycles per second. This may be the reason why 
plots of log 7 vs. 1/T show anomalous behavior 


7. 


which Stambaugh interpreted as due to crystal- 
lization phenomena. 


Breaking of Polymeric Threads Under Load 


In certain recent experiments on cotton and 
rayon threads,”° the life of these threads at 
various temperatures was measured as a function 
of the load. 

Suppose that the breaking is due to the slipping 
of bonds (which act as “primary” bonds for 
small loads, but which may be micellar bonds). 
If the number of such bonds per unit area is N, 
the rate of breaking of such bonds under stress fp 
if repair is possible is 

1dN kT 


Ndt kh 
fa 
Xexp (—AF*/RT)2 sinh 
2Nk71 





(25) 


For large values of the stress, or when repair 
of the bonds does not take place, the equation is 


1 dN kT 


Ndt kh 
Xexp (—AF*/RT) exp (fodA/2NRT). (26) 


0 Busse, Lessig, Loughborough, and Larrick, J. App. 
Phys. 12, 715 (1942). 
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TABLE I, 
. dynes 
- Polymer n3 poises GitGs cm? G2/r2 
Natural gum 3 X10 2.8X 107 3.67 X 10° 
Natural rubber 1.4 104 9.6 107 4.2 x10? 
Synthetic D 3.5 X 104 1.7 108 1.63 X 108 
Thiokol 3.2 104 1.4 108 1.0 108 
Neoprene GN 2.6 10! 1.4X 108 7.5 107 
3 large values of x 
5 e 7 
S — Ei(—x)2—. (28) 
x x 
g 
g Equation (27) therefore becomes approxi- 
3 mately 
c 
= e-fO2NKT pT 
a ————— =—- exp (—AFF¥/kT)r. (29) 
for/2NokT h 
Or rewriting 
— for for kT AFF 
log r= —log —log —+——. (30) 
4 2NokT 2NokT h RT 





a 
° i i i 
° 5 to) 1S 20 25 30 


1/v? X 104 (C.P.S. at 28°C) 


Fic. 8. A, rubber gum stock; X, rubber; OH, synthetic D; 
o, thiokol; ©, neoprene E and GN. 








This equation is in essential agreement with 
the experimentally observed equation of Busse, 
Lessig, Loughborough, and Larrick. It is to be 
noted that: 


‘ . 2 d log T for Navogadro 
By the substitution y=foA\/2NRT we obtain R——— = AH* ———_—__-. (31) 
the following: Pid 2No 
. - L 
Bat =~ exp (—AFF/kT)t, (27) It is to be-noted that [Rd log r/d(1/T)] is 


where 7 is the lifetime of the thread and Np» the 
initial number of bonds per cross section. The 
upper limit of the integral is infinite because 
when the thread breaks the number of remaining 
bonds is zero. 

The integral is the well-known and tabulated 
logarithmic integral —Ei(—x). For sufficiently 





only the apparent heat of activation. To obtain 
the true heat of activation one must either 
extrapolate the apparent heat to zero load or add 
(forANavogadro)/2No. 

From the data of the above paper it appears 
that N, the number of bonds per unit area, is 
about 7X10", and AFF, the free energy of 
activation, is about 27 kcals. 
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The Diffusion of Air Through Monel Metal 
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The rate at which gases from the atmosphere diffuse through monel metal at various temper- 
atures has been observed. An evacuated heavy-walled metal tube, a section of which was 
machined to a thin wall, was heated in a furnace. The rate at which the pressure inside the 
tube increased as a result of the diffusion of gases from the atmosphere was noted. The log R 
vs. 1/T curves, which are usually linear, show breaks about 900°C indicating some unusual 
physical change occurring in the alloy at this temperature. Dilatometric measurements on 
the same alloy were used to check the location of this break. 


REVIOUS studies of the diffusion of gases 

through metals have been largely confined 
to hydrogen-pure metal systems. A few studies 
of the diffusion of oxygen and nitrogen through 
certain metals have been made. These studies 
have been carefully summarized by McBain,! 
Smithells,? and Barrer*® and, therefore, will not 
be reviewed here. 

It is the purpose of this paper to present data 
relative to the diffusion of gases from the 
atmosphere through a common alloy of nickel, 
copper, and iron. This alloy is of great importance 
because of its wise usage in research apparatus, 
its non-corrosive properties, and the ease with 
which it may be machined at high speeds. 

Briefly the experiment to be described is: 

A metal tube closed at one end, the other end 
being connected to a vacuum system was heated 
in a furnace. The diffusion of gases from the 
atmosphere was observed by noting rate at 
which the pressure increased on the inside of the 
tube. The diffusion rate was observed at various 
temperatures. 


THE APPARATUS 


The apparatus essentially consists of the 
diffusion tube, the furnace, a vacuum system, 
and a temperature control. These parts are 
shown in Fig. 1. 


The Diffusion Tube 
The diffusion tube was made from heavy- 


1J. W. McBain, The Sorption of Gases by Solids (Rout- 
ledge, 1932). 

2C. J. Smithells, Gases and Metals (Wiley, 1937). 

3R. M. Barrer, Diffusion in and Through Solids (Mac- 
millan, 1941). 


walled tubing (4.5 mm), a section of which was 
machined to a wall thickness of 1.15 mm. The 
inside diameter of the tube was 1.74 cm. The 
length of this section of tube was 10 cm. One 
end of the tube was closed by welding a nickel 
disk on to the tube. Atomic hydrogen welding 
equipment was used. The tube was obtained 
from The International Nickel Company, and 
was found to be of the following chemical 
composition: 69.70 percent nickel, 28.00 percent 
copper, and 2.57 percent iron. 

In addition, carbon was found to be present 
and estimated to be 0.1 percent. 


The Furnace 


A Nichrome resistance furnace was used to 
heat the diffusion tube. This furnace was 
provided with special booster coils near its 
open ends to reduce temperature gradients. 
This is shown in Fig. 2. 


The Vacuum System 


The open end of the diffusion tube was 
connected to a vacuum system by means of a 
glass to metal joint using DeKhotinsky cement. 
This joint was kept cold by means of a water 
cooler as indicated in Fig. 2. 

The vacuum system was arranged with suit- 
able stopcocks such that the diffusion tube and 
McLeod gauge could be isolated. A mercury 
diffusion pump and a Cenco Hi-Vacuum pump 
were used to exhaust the system. See Fig. 3. 
The fore pump was arranged with a special 
exhaust which made it possible to collect the 
diffusing gas by the displacement of water. 
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The Temperature Control 


The temperature was measured by means of 
two Chromel-Alumel thermocouples fastened to 
the outside of the diffusion tube. One thermo- 
couple was placed near the center of the thin 
section of the tube, the other at the edge as 
indicated by Fig. 2. The relative currents in the 
main heating coil and the booster coils were 
adjusted to reduce the temperature gradient 
along the length of the thin section of the 
diffusion tube to a minimum. The temperature 
was controlled both manually and automatically 
using a capacitrol in conjunction with a Leeds 
and Northrup portable potentiometer. See Fig. 4. 
The temperatures did not deviate from the mean 
by more than 2°C. 


PROCEDURE 


After the diffusion tube was evacuated and 
known to be vacuum tight, it was heated to 
1250°C and degassed for one hour. Then a series 
of isothermal runs were made by closing the 
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stopcock leading to the pumps and noting the 
rate at which the pressure built up inside the 
diffusion tube as a result of the diffusion of gases 
from the atmosphere through the thin section of 
the diffusion tube. The stopcock leading to the 
pumps was opened between each run to remove 
the diffused gas at which time the temperature 
was changed for the subsequent run. This 
technique of observing diffusion is known as 
the “trapping method’ and is described in 
detail by Ham and Post‘ in their report on 
hydrogen diffusion through nickel and iron. Also 
the “‘steady state’’ method of observing diffusion 
developed by Ham*® was used to check the 
trapping rates. The ‘‘steady state’’ method 
involves the partial closing of the cock leading 
from the diffusor to the pumps. The cock is 
adjusted by partly closing it until a desired 
pressure builds up in the diffusion tube. This 
pressure is permitted to come to equilibrium so 
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4C. B. Post and W. R. Ham, J. Chem. Phys. 6, 598 
(1938). 
5W.R. Ham, Trans. A. S. M. E. (June, 1937). 
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that the pumps are removing the gas at the 
same rate at which it enters. The magnitude of 
this build-up pressure depends upon the temper- 
ature and has been shown to be proportional to 
the actual diffusion rate by Ham and Post.‘ 
The results of this experiment tend to confirm 
their claims for the ‘“‘steady state’’ method. 


THE DATA 


The data are presented in graphic form using 
the usual coordinates of logio R vs. 1/T, where 
R is the diffusion rate and T is the absolute 
temperature at which the diffusion rate is 
determined. 

(1) Figure 5 shows the degassing rate of the 
tube using the “‘steady state’’ method during a 
temperature range from 20°C to 1250°C and 
from 1250°C to 400°C. 

(2) Figure 6 shows the diffusion rate as a 
function of temperature. The temperature was 
varied in a cycle from 450°C to 1250°C and from 
1250°C to 450°C. The “‘steady state’? method 
was used in this run. 

(3) Figure 7 shows the trapping curves for a 
temperature range of 672°C to 1150°C. The 
pressure is expressed in microns of mercury and 
the time in minutes. 

(4) Figure 8 shows the relation of the pressure, 
built up during the first minute after the gas 
was trapped, as a function of temperature. 

(5) Figure 9 shows the steady state data which 
were obtained as a by-product in obtaining the 
trapping data. Each time the temperature was 
increased for the next run, the cock which iso- 
lated the diffusor from the pumps had to be 
opened in order to remove the diffused gas. 
Even with the stopcock -wide open, there is a 
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slight pressure built up in the system. This 
“steady state’”’ pressure is a function of temper- 
ature and is referred to as the zero sensitivity 
steady state rate. 

(6) Figure 10 shows the functional relation 
between the pressure and time for the trapping 
data. The plots are made using the coordinates 
P* vs. t. The pressure P is expressed in microns, 
and the time ¢ in minutes. This method of 
plotting is used to rectify the P—¢ curves and 
make possible the determination of the velocity 
constants at various temperatures. 

(7) Figure 11 shows a logarithmic plot of the 
slopes of the rectified pressure-time curves as a 
function of the reciprocal of the absolute 
temperature. 

(8) Figure 12 shows “steady state’ data for an 
entirely new diffusor and completely rebuilt 
apparatus used to check the results indicated in 
Figs. 5-11 inclusive. 
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(3) The slopes of the log R vs. 1/T are variable 
in certain regions. These regions are referred to 

(1) It is evident that gases from the atmos- as ‘“‘breaks.”’ 
phere diffuse through an alloy of nickel, copper, (4) The breaks in the log R vs. 1/T curves for 
iron, and carbon. the various runs occur in the same temperature 

(2) The slopes of the log R vs. 1/T curves region. 
obtained by the “steady state’’ and by the (5) The exact cause of the breaks in the 
‘‘trapping methods” are substantially the same. log R vs. 1/T curves is not known. It is possible 
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that it is a result of either the formation of a 
solid solution of NiO and Ni or an order-disorder 
occurring among the constituents. 

(6) The breaks in the log R vs. 1/T curves are 
reproducible and do not appear to be greatly 
influenced by the heat treatment received during 
the experiment. 

(7) The break in the Jog R vs. 1/T curves 
shows a hysteresis very similar to that obtained 
by Ham and Post! in their study of hydrogen 
diffusion through nickel and iron in the Cure 
regions for these metals. 

(8) The material of diffusion tube was tested 
magnetically and was found to be non-magnetic 
above 50°C. Therefore the break occurring at 
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900°C in the log R vs. 1/T curves could not have 
been a magnetic change. 

(9) The functional relation between the pres- 
sure and time of the reaction rate curves indi- 
cates that the diffusing gas leaves the inside 
surface of the tube as a diatomic molecule. 

(10) The mechanism by which gases from the 
atmosphere diffuse through the tube is not known 
except that the gas is an atomic form as indicated 
by the reaction rate curves. 

(11) The break in the log R vs. 1/T curves is 
an indication that there is some physical change 
occurring in the alloy at about 900°C. This 
results in a change in resistance of the alloy to 
the passage of the ironized atmospheric gas (or 
gases), much like an electrical resistance change 
to an ordinary electric current. This being the 
case, one would expect to find deviations in 
other physical properties at 900°C, such as 
thermal expansion, density, etc. Thermal ex- 
pansion data were obtained by dilatometric 
methods using a section of the actual diffusor as 
a sample. These data are shown on Fig. 13 and 
indicate that there is a sharp change in the 
coefficient of expansion at about 925°C. This is 
an excellent check on the break obtained by 
diffusion methods. We are indebted to Dr. Carl 
H. Samans of The American Optical Company 
for the dilatometer curves. 
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A Note Concerning /-Type Doubling in Linear Polyatomic Molecules 


HARALD H. NIELSEN AND WAVE H. SHAFFER 
Mendenhall Laboratory of Physics, Ohio State University, Columbus, Ohio 


(Received January 4, 1943) 


It is shown that certain Coriolis interaction terms in the rotation-vibration Hamiltonian 
of a linear polyatomic molecule give second-order contributions to the energy which remove 
the degeneracy of levels associated with the quantum number / of internal angular momentum 
arising from a doubly degenerate perpendicular mode of oscillation. This type of splitting 
is called “‘/-type doubling” because it is exactly like ‘“‘A-type doubling” in the electronic energy 


levels of diatomic molecules. 





TYPE of oscillation occurring in linear 

polyatomic molecules is one in which the 
constituent atomic nuclei are displaced perpen- 
dicularly to the internuclear axis of the molecule. 
Such oscillations are isotropic in two dimensions 
and are known as perpendicular modes. The 
quantum states of such oscillations are defined by 
two quantum numbers »; and J,; the first of these 
is the total vibration number which may assume 
positive integral values and the second is a 
measure, in units of /27, of the internal angular 
momentum directed along the internuclear axis 
which is associated with the doubly degenerate 
mode of oscillation. For a given value of 1; the 
quantum number J, may assume the values 
+uy,, +(%,—2), +(v,—4), ---+1 or 0. In zero 
order of approximation the vibration energy 
E,.°=hcw,(v~z+1) is independent of J,; hence such 
a state is (v,+1)-fold degenerate in 1];,. The 
rotation energy depends on /,? and in seeond 
order of approximation the anharmonic terms 
in the potential function of the molecule also 
contribute to the energy terms which depend on 
l.2. The degeneracy in 1, is, therefore, partially 
removed up to second order of approximation 
and the zero-order vibration energy level splits 
up into (v,+1)/2 components for odd values of 
v, and into 1+ (v,/2) components for even values 
of v,. The energy remains independent, however, 
of the algebraic sign of ), and each sub-level 
remains twofold degenerate. 

In a recent paper Herzberg! has suggested that 
there must be terms of interaction between the 
rotational and vibrational motion present in the 
Hamiltonian for a linear molecule which will 
further remove this degeneracy. In an attempt 


1 Gerhard K. Herzberg, Rev. Mod. Phys. 14, 219 (1942). 


to clarify the character of these perturbations he 
has presented certain arguments of a semi- 
classical nature. In particular he has pointed 
out that certain terms in the Hamiltonian for the 
linear X2Y2 type molecule developed by Shaffer 
and Nielsen? (SN) have matrix elements diagonal 
in the quantum number y, but non-diagonal in J,. 
These terms lead to a splitting of the levels and 
Herzberg has computed the magnitude of this 
splitting for T[-vibration states (i.e., where 
},=1). Herzberg has, moreover, recognized the 
similarity of this type of splitting to that known 
as A-type splitting in diatomic molecules and he 
has accordingly termed it ‘‘J-type doubling.’”’ He 
has, moreover, demonstrated that such an /-type 
doubling is consistent with the existing data on 
the spectra of CO2 and C:He. 

While we cannot subscribe to the classical 
arguments of Herzberg except in a very rough 
way, it is indeed true that there exist interaction 
terms in the Hamiltonian for a linear polyatomic 
molecule which further remove the degeneracy 
in J,. The purpose of this note is essentially to 
point out that the Hamiltonian for the X2Y:2 
linear molecule developed by SN, although 
entirely adequate to give the elements of the 
energy matrix diagonal in v and / to second 
order of approximation, is not suitable for com- 
puting the correct value of the /-type splitting. 
This is not due to an error in computation, but 
as the authors indicate, certain terms which did 
not contribute to the diagonal elements of the 
energy matrix in this approximation were 
omitted to save space. We simply state here the 
correct form of the Hamiltonian for this type of 


2W. H. Shaffer and A. H. Nielsen, J. Chem. Phys. 9, 
847 (1941). 
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molecule and then compute the value of the X2Y»2 molecule has the following form complete 
l-type splitting. The calculations could be carried to second order of approximation and suitable 
out readily for the general linear polyatomic for calculation of the /-type splitting: 

molecule; we have chosen the linear X2Y2 model . 

because of its relative simplicity and especially j7o= (he/2){¥ wi(p2/h?+q2) 


since it is one of the models considered by i=1 
Herzberg in his recent paper. . 2 oe 2 2 
After the usual contact transformation has +h wnL (Dirt pie) /h?+ (quit dee) J} 
been performed the Hamiltonian for the linear , +(1/2I,)(P2+P,?). (1a) 





H® = (h?/2I,) —(th/Ie)_8°qipitv7q2h2— v5(q2p1 + Gip2) ]+ (1/21.) (p22 + py) + U2 
— (1/2) (Ze)8L (8?/A1) + (v?/A2) (P22 +P,*)?+ (1/20) (aiigi? +411’ pi?/h? +2092" 
+022! p2?/h? +3393" +33’ ps*/h? +44( Gar? +G42”) +44! (Par? + Pas?) /h? +055(G51? +952”) 
ass (Psi? + ps2?) /h?} (P2+P,?) + (1/Te) {baa(qar’P2+garP,?) +44 (pas?P2+ par*Py?) /h? 
+b55(qs2?P 2? +951°Py*) +55’ (ps2*P 2? + psi2P,?)/h?} — (1/T-) {basqarga2 t+ baa Parpa2/h? 
+ 55951952 +555 PsiPs2/h?} (P2:Py+P,P.) —he} z (B:i:?/w,)(1+(39:4/2) + (3p,79:7/h?*) J 


2 3 
+0 DL (Bij;?/20:) (40? —w?) [| (20? — @,*)g;' +400 59.7q7? + 207977 7/h*] 
i=1 i~j=1 


2 5 
+>. > (Bin, 2w ;) (404? — wi?) "[ (Qe? — wi?) (Gur? +2”)? +40 ong i? (Gur? + x2”) 


i=1 k=4 


2 3 2 5 
+ Qe? (Gir? Pir? + Guo” Pu2”) /h? | +3 } ps (BiB ij;/wiQi7q7+3 pm >» (BiB ins/wigi? 


i=1 j=l xi i=l k=4 


2 3 5 
X (Ger? + x2”) + (B1228133/w1)q2"GQs? + (82118 233/w2)gi7gr? +>, DY YX (BijiBinx/wiq? 


i=1 i4j=1 k=4 


X (ger? 922) +L (B i448 is5/wi) (G41? +42?) (Go1? +952”) }. (1b) 
i=1 


The notation used here is that of SN except for the substitution of gi: for r, and qx2 for s,. The fol- 
lowing constants are defined differently from the corresponding ones in the paper of SN: 


diy = (A/T dy) {3 —2y?(A1 — 2X5) /(A1— As) } 
+4mc(h/I.)*\ (281115/d1*) — Birzy(3A1 — Az) /[A2? (441-2) J}, 
yy’ = 2(h/Ted1*)y7Ai/ (Ai — As) +4-c(h/Te)* | (B1118/A1*) — (BrrzyA1/[A2?(4A1—Az) J}, 
d22= (h/ Ido") {3 —262(A2— 2As)/(A2—As) } 
+4mc(h/I.)'{ (—2B2227/A2*) +B81226(3A2—A1) /[A1*(4A2—A1) J}, 
doo = 2(h/ edo!) 5°A2/(A2—As) +4ac(h/T.)§{ (—Bovxy/A2*) +B1225d2/[A13(4A2—A1) J}, 
33 = (/ Teds) (As+Aa)/(As—Ya) 
+4arc(h/I.)*{ 6\1-*B331(3A3 —A1) / (4X3 —A1) — YA228332(3A3—Az) /(4A3—Az)}, 
33’ = 2(h/Ids*)A3/(As— Aa) +4arc(h/T)*{ 5A1-*Bg31A3/(4A3—A1) — YA2~*B332A3/(4A3—Az2)}, (2) 
44 =4rc(h/Te)*{ 5B aard1~2(3X4—A1)/ (44 — As) — VBaa2do72 (3X4 — Az) /(4A4—Az) }, 
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O44’ =4mc(h/Te)*{ 5Bast1*Aa/(4A4—A1) — VBaa2d2*Aa/(4A4—Az)}, 
A55=4mc(h/I.)*{ 5Bs51A1~* (3X5 —A1)/(4A5s—A1) — VBss2d27?(3A5—Az)/(4A5—Az)}, 
55 =4c(h/I.)*{ 5Bss1A1*As/(4A5—A1) — VBss2d2*As/(4A5— Az) }, 

baa = (h/2T.)rs*(Aa+A3)/(As—A3), 
baa’ = (h/Te)da*da/ (Xa — Hz), 

bss = (h/2T.)As*{ (As +A2)/(As—Az) J6°+[ (As +A1)/(As—Az) }77}, 
bss’ = (h/2T.)s*{ [Ns/(As— Az) 5? J +[A5/(As—A1) Jy}. 


As Herzberg has pointed out, it is just the terms containing @44, @44’, Das, baa’, G55, G55’, etc., which 
will give rise to the /-type splitting. These terms originate entirely with the Coriolis forces and are 
not, as Herzberg suggests,’ in part due to the harmonic deformation of the molecule. When one 
replaces gi by px COS dx, Giz by px Sin dx, tPkr/h by 


fe] sin od: 0 
( cos ¢,—— —) ; 
Opr pr Ody 





and 1px2/h by 











} 0 cos¢ 0 
(sin oi— + —) . 
Op, = pe OH 
the term in question reduces to 
5 0? 1 ¢@ 1 3 
(1/2J,.) i busn?—bur| ~~ —mepen |[eetra 
k=4 Op.” px Op, px” Ady?” 





- a2 10 1 @ 
— (1/41) & ef bape? bul] —+— — + — | 
k=4 
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Pr Op.Od, pr? Ode =k Ope =p.” Oy” 
—_— a2 108 1 @ 
—(1/4I.) 2 e+ bupst—bu| —+— ———— | 
bane Op.” pr Op; px” Ady” 
yl #@ 10 t 0 ¢ @ : 
+2ibw'|— He — ltr Py+iP.P+ P,P.) (3) 
Pr Opr.Ode pr? Ode pK Apne MK” AG,” J 








lowing: 


(val 4vsl 5 | H® |v4' 1 4'05 Ls’) 


The matrix elements of H diagonal in all the 
quantum numbers, including the first term in 
(3), have been calculated in the work of SN and 
will not be repeated here. The additional matrix 


— ¥ oy * O51 J QDYog v "- 
elements which we shall investigate are the fol- f Vous eos R* HO YeuousRsvdr, (4) 


_ 3That these terms originate entirely with the Coriolis where v,v,’ and 1,41,’ and where 
forces may be seen from an inspection of a problem sug- 


gested to us by Herzberg in a private communication. The 











quantum mechanical Hamiltonian for the linear XY2 
molecule is derived from a classical kinetic energy expres- 
sion where the classical Coriolis terms have been set equal 
to zero. The eigenvalues are readily obtained and one sees 
that the J-type doubling is identically zero. 





Porn = Nox exp ((— pe?/2) +41 cy) Fs (on +n) (pn?) 
if | 
Nowe = 24{[ (ve —1e)/2]'}4/ (Le + 1a) /2 443, 


t) 








and F*s(.+) is an associated Laguerre poly- 
nomial, and R,; is the rotator wave function in 
which 1/= >>; I. 

To calculate the energies arising from /-type 
doubling it is necessary to evaluate only those 
matrix elements for which v=,’ and 1,#];,’ 
because any element not diagonal in v, can con- 
tribute to the energy only in a still higher order 
of approximation. The second and third parts of 
(3) are the only terms in H® which have non- 
vanishing matrix elements when »%,=2' and 
lL. 1,’, and only for 1,’ = l,+2, because the only 
non-vanishing elements of P,’, P,?, P,P,, or 
P,P, are those for which either ],=/,’ or 
L,.’=h,.42. It can be shown that 


(vex | px? | Vel, 2) = C (vl) (v,-+1.+2) ]} 


and that the corresponding matrix elements of 
0? 1 ¢a 1 @ ) 
Op.” pr Opx px” Ody.” 


have the same values as those of p,2. It can be 
shown that 





Pr Oped, px? Oe pe Apx px” Ox” 





(- 0 1d +2 6+ # 


has no non-vanishing matrix elements diagonal 
in v. Integration over the rotational wave func- 
tions gives 


(J1| P,2|Ji+2)=—(J1|P,?| Jl +2), 
i 
= +-(J1 |P,P,+P,P,|J1+2), 


= (h?/4)(f—1(1 +1) PLf—(141)(1 +2) }', 
where 
f=J(J+1) and l=>.k. 


The non-diagonal matrix elements of H® which 
contribute to the /-type doubling have the value 


(y%1,| 1 | od. 2) 
= (hcB2/on){ (Oe i) (He tMe +2) }} 
XCf—1(d +1) PLf—(241)(1+2)]', (5) 


where B,=(h/8°cI.). The values of the energy 
associated with the /-type doubling can be deter- 
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mined from the secular determinant 
| (vl), | TT® | ode’) ae Eby. | =(), 


The twofold degeneracy in /;, associated with the 
independence of energy of algebraic sign of I, 
is thus removed not only for /,;=1, as Herzberg 
indicates, but for all values of /,. A few examples 
of l-type doubling are considered below. 

For v,=1 the quantum number J, takes the 
values +1 and —1 and if the quantum number 
l,»=0 we have what may be termed a II-vibra- 
tion state. In this case the only non-vanishing 
off-diagonal matrix elements are 


(1, +1!) H@|1, 1) =2hce(B.?/o,)J(J+1). (6) 
The two roots of the secular determinant are 
E=E+2hc(B,2/wx)J (J +1), (7) 


where E denotes the energy value given by SN. 

When »;,=2 the quantum number J, may take 
the values +2 and 0, and if l,.-=0 we have 
what may be termed a A-vibration state and a 
y-state. The non-vanishing off-diagonal matrix 
elements are 


(2, +2|H®|2, 0) =(2, 0|H®|2, +2) 
= 2he(B.*/wx) f(f—2) =A. 


If EZ, and Ey denote for ,=+2 and i,=0, 
respectively, the values of the energies given by 
SN, the secular determinant for the /-type 
splitting for this case is 


E.-E A 0 
A E,-E A |=0 (8) 
0 A E.-E 


The roots are E= EF, and 
E=(1/2){(2o+E2)+[(B2—E)*+8A*]}!}. (9) 


By referring to the paper by SN the reader may 
readily verify that (E,—E,)=4(Xiu—B.). The 
magnitude of (E.—E,) is, in general, large com- 
pared with (B,?/w,;)? and the latter two roots in 
Eq. (9) may be expanded and expressed approxi- 
mately as 


E=E.4+f(f—2)(4h?c?B.'/w,2)/(Xizn—B.), (10a) 
E=E,y—f(f—2)(4h?c?B,4/w,?)/(Xix—B.). (10b) 
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The /-type doubling thus produces a splitting of 
the A-state of magnitude (4hcB,4/w,”) f(f—2)/ 
(Xix—B.) cm, that is, the splitting is propor- 
tional to (B.?/w,)?/(Xix—B-). Ina similar manner 
it may be shown the /-type splitting for any 
vibration state is proportional to (B.?/w,)' 
X (Xi. — B.)~"**1. It may be verified by reference 
to the works of Hill and Van Vleck‘ and Kronig® 
that the /-type splitting in the case of funda- 





mental and harmonic vibration states is entirely 
analogous to A-type splitting for singlet elec- 
tronic states in the case of diatomic molecules. 

An interesting case of /-type doubling may be 
found in the example of the combination state 
wa+ws. In this case v45=v5;=1 and J, and /; may 
each take the values +1. If a=2hc(B,?/w,)f, 
b=2hc(B2/ws)f and f=J(J+1) the secular 
determinant for this case is: 

















I,=i,=1 —l4=1;=1 —ls;=1,=1 l4=1l,=-1 
1,=l,=1 B.—E a | b | , 0 
= sul =1 a E.-E 0 | b 
. | =0, (11) 
—l,;=14=1 b 0 E,-E | a 
l4=ls=-1| 0 b a | B,-E | 














where FE, and E, are obtained from the energy 
expressions of SN by setting /4=/;=+1 and 
l4=—ls;=+1, respectively. The four roots of 
Eq. (11) can readily be shown to be, since 
E.—E,)=4B,, 
E=E.+hc?B,?f2(wstzws!)?, (12a) 
E=Eo+h?c?B.3f?(ws ws 2. (12b) 


The /-type splitting for this case thus becomes 
4hc(B.°/wws)f? cm. This case resembles the 
A-type doubling of doublet states in Hund’s case 
a. 

The preceding considerations may as readily 
be applied to the general linear polyatomic 
molecule discussed by Nielsen.’ In addition to 
the matrix elements diagonal in all the quantum 
numbers given in that work it is necessary to 
consider the off-diagonal elements corresponding 
to those given in Eq. (6) for the linear X2Y2 
model: 

(vi1;- ° “Vel xe ‘Unb n| H® loyly- ° Vel, 42° . *Unl n) 
= D> ofes2hc(B.?/ wx) (uF Ly.) }(Op L,.+2)3 
X(f—1(0 +1) PLf—(7£1)(7 +2) }', (13) 
4E. Hill and J. H. Van Vleck, Phys. Rev. 32, 250 (1928). 


5 R. de L. Kronig, Zeits. f. Physik 50, 347 (1928). 
6H. H. Nielsen, Phys. Rev. 60, 794 (1941). 











where / = >>; 1. The expression (13) differs from 
(5) only in the presence of the coefficient >> sé3? 
where &s is one of the Coriolis parameters and 
depends upon the masses of the nuclei, the dimen- 
sions of the molecule and the harmonic force 
constants. In all the linear cases where we have 
evaluated the sum >>; &? it has been equal to 
unity. This is probably true for the general linear 
molecule but we have no general proof of it. 

While the /-type doubling exists for all states 
having a quantum number /;, 0, it is likely that 
only in II-vibration states is the effect large 
enough to be of importance in the interpretation 
of spectra. 

It may be of interest to remark that /-type 
splitting must also occur in the perpendicular 
vibration states of axially symmetric molecules. 
This problem is very complicated and not worth 
considering since bands of this type are in general 
so complex that the effect could not be observed. 
Since the linear molecule is a special case of the 
axially symmetric molecule, the problem of /-type 
splitting in axially symmetric molecules could be 
handled by a procedure similar to that given in 
this paper. 
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Crystal Structure of Gadolinium Formate, Gd(OOCH); 


A. PABST 
University of California, Berkeley, California 
(Received December 8, 1942) 


The crystal structure of gadolinium formate has been studied by means of Laue, rotation, 
and powder x-ray patterns. It has a rhombohedral lattice. The unit cell, a9 =6.17A, a=115° 30’, 
contains one molecule, Gd(OOCH);. Density—(calculated) 3.85, (Berman balance) 3.77. 
By taking into consideration the symmetry of special positions all space groups except D;7— R32 
and C3,5—R3m can be excluded. A structure has been found in the latter group which yields 
excellent agreement of calculated and observed intensities. 


\v 


\v 





INTRODUCTION AND ACKNOWLEDGMENTS 


EARLY ten years ago the writer reported 
observations on crystals of gadolinium for- 
mate.' Sometime later more crystals of this sub- 
stance were kindly prepared by Dr. George C. 
Nutting, then of the College of Chemistry at the 
University of California, and these have been 
used for the determination of the crystal struc- 
ture. Most of the x-ray diffraction patterns used 
in this work were obtained while the writer held 
a fellowship of the John Simon Guggenheim 


TRANSLATION LATTICE 


Powder patterns showed only lines for which 
h—l—k=3m, indicating a rhombohedral lattice. 
The rhombohedral unit cell, containing just one 
Gd(OOCH);, has the dimensions a,=6.17A, 
Za=115° 30’. 

The morphology of the crystals’ is clearly indi- 
cative of the lattice character. The following 
table shows the agreement with the theoretical 
order of forms.‘ 


, a ; ‘ Forms Order 
Memorial Foundation in the Department of : 
: os # {hkil} {h,k,l,} Observed Donnay-Harker law 
Mineralogy, British Museum of Natural History. — ; 4439) (101) ; \ 
The kindness of Dr. W. Campbell Smith and of {1071} {100} 2 2 
Mr. F. A. Bannister in making available the {1010} {112} — 3 
{0221} {111} 3 4 


facilities of that department is gratefully ac- 
knowledged. 


DIMENSIONS, DENSITY, AND CELL CONTENT 


From a powder pattern with Mo radiation in an 
8-inch cassette, later checked by numerous other 
patterns, the dimensions of the smallest possible 
hexagonal cell were found to be a,=10.44A, 
Cr=3.98A, c,/an=0.381. The axial ratio pre- 
viously reported,’ 0.375, is probably less accu- 
rate. It was based on the angle p of the unit 
form, {1011}. The axial ratio corresponding to 
the observed angle p of the {0221} form is 0.381, 
just checking the x-ray value. 

The density of gadolinium formate crystals 
was determined by Berman torsion microbalance 
to be 3.77+0.04. This agrees satisfactorily with 
the value 3.85, which would correspond to a 
content of three formular units in the smallest 
hexagonal cell. 


1A. Pabst, Am. J. Sci. 36, 72-79 (1933). 
* Reference 1, p. 73. 


SYMMETRY AND SPACE GROUP 


A Laue pattern with the beam normal to the 
(101) face shows that the crystal belongs in the 
Laue group D3a—3m. The only extinctions are 
those required by a rhombohedral lattice. Rota- 
tion patterns with Cu and with Cr radiation 
show every line, even including all beta-lines, 
permitted by a simple rhombohedral lattice. 
This indicates the space groups D3,.°—R3m, 
D;'— R32 or C3,5>—R3m. 

The choice of space groups can be further 
limited by considering the symmetry of special 
positions. The rhombohedral unit cell contains 
three C atoms. From the known configuration of 
the formate group these C atoms cannot lie at 
centers of symmetry. In the space group R3m 
there are two sets of threefold positions, both at 

3 Reference 1, p. 72. 


4]. D. H. Donnay and David Harker, Am. Mineral. 
22, 446-467 (Palache number) (1937). 
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centers of symmetry. This leaves no suitable 
position for C and eliminates the group. In the 
groups R32 and R3m the threefold positions are 
on twofold axes and planes of symmetry, respec- 
tively, and either kind of position is compatible 
with the configuration of the formate group. 

Auxiliary observations sometimes help in the 
choice of space groups which lie in different sym- 
metry classes : 

Morphology in this case® suggests the class 3m, 
just eliminated, so that the crystals are certainly 
pseudosymmetrical. A statistical study of re- 
corded faces on eleven’ crystals gives no indica- 
tion of a tendency to hemimorphism which 
might be expected in class 3m. 

Pyroelectricity is a symmetry criterion in this 
case only if the distribution of charges on the 
surface can be observed with exactness, since 
both possible classes, 32 and 3m, would permit 
the effect. The small size of the crystals, diameter 
less than one millimeter, precluded the possi- 
bility of satisfactory observations. 

Optical Activity would indicate the class 32. 
It was not detected, but this does not eliminate 
the class since the effect generally is very slight 
and not easily observed. 

Etch Figures developed on {1120} {101} indi- 
cate no symmetry normal to that form. This, 
taken together with the morphology, was pre- 
viously interpreted® as indicating the class 
C3;—3, but that class is incompatible with the 
Laue symmetry or the requirements of the car- 
bon positions. In class 32 these figures should 
show twofold axes. Since they do not they might 
be considered to eliminate it. Because the etch 
figures are not very good it seems safer to say 
that they merely yield a presumption in favor 
of class 3m. 


DETERMINATION OF THE STRUCTURE 
The positions in space group R32 are: 


1: (a) 000 (b) 334 

2: (c) xxx; £2% 

3: (d) Oxt; 20x; x%0 (e) $x; Zhu; xE4 

6: (f) xyz; sxy; yox; G22; 292; B27. 
The choice between positions (a) and (b) corre- 
sponds merely to the selection of an origin. It is 


5 Reference 1, p. 72. 
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required to assign coordinates to: one gadolinium, 
three carbon, three hydrogen, and six oxygen 
atoms. Assigning Gd to position (a), the carbons 
may be placed on the twofold axes in (d) or 
in (e). It will, of course, be impossible to deter- 
mine parameters for hydrogen from intensities, 
but the hydrogens will lie on the same axes “is 
the carbons. All the oxygens may then be placed 
in a sixfold ‘“‘general’”’ position. This allows the 
placing of six oxygens about each gadolinium 
in a slightly deformed octahedral arrangement, 
each formate group being tied to two gadolini- 
ums. Several variants of this structure may be 
found which conform to the known configuration 
of the formate group® and the radius of gado- 
linium. A certain variation of the parameters is 
possible, sliding the formate groups and rotating 
them on the twofold axes, within the limits al- 
lowed by atomic radii. Intensity calculations 
were carried out for numerous combinations of 
parameters within these allowable limits. None 
of them showed satisfactory agreement with ob- 
served intensities and any structure in this space 
group must be considered eliminated. 
The positions in space group R3m are: 


1: (a) sxx 
3: (b) xxz; x2x; 2xx 
6: (c) xyZ; 2xy; yox; yxs;.syx; xZy. 


The choice of a parameter for gadolinium in the 
(a) position corresponds to selection of an origin 
and it is best placed at 000. The carbon atoms 
necessarily go into (b) positions which lie in the 
planes of symmetry. If all the oxygen atoms are 
equivalent they must lie in (c) positions, flanking 
the symmetry planes and forming ditrigonal 
rings, which might be practically hexagonal, 
about the threefold axes on which lie the gado- 
liniums. Such an arrangement does not permit 
the tying together of the structure by formate 
bridges between gadoliniums because the two 
oxygens of each formate group would both be 
closest to the same gadolinium, as may be seen 
from an inspection of the symmetry diagram for 
this space-group.’? Such a structure would be 
molecular and possibly even sheetlike. From the 


6 L. Pauling and L. O. Brockway, Proc. Nat. Acad. Sci. 
20, 336-340 (1934), and W. H. Zachariasen, J. Am. Chem. 
Soc. 62, 1011-1013 (1940). 

7 International Tables for the Determination of Crystal 
Structures (New York, 1935), vol. I, p. 243. 
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_Fic. 1. Top: Section through rhombohedral cell of gadolinium formate in the 
1210 or 011 plane, which is a plane of symmetry. The Gd atoms lie on vertical 
threefold axes which are the intersections of three such planes. Scale is given by 
length of rhombohedral axis indicated in diagram. Parameters as at head of 
column 5, Table I. Bottom: Sections through half cells, one-half scale of top 
diagram. Showing, from left to right, the orientation of formate groups corre- 
sponding to parameters in columns 3, 4, and 5 of Table I. 


properties of the crystals, hardness, brittleness, 
prismatic habit, and absence of cleavage, this is 
clearly not the case. 

The final alternative is to place the formate 
groups in the planes of symmetry, the oxygen 
atoms going into two non-equivalent sets of 
threefold positions. The range of positions com- 
patible with atomic radii and the structure of 
the formate group can be easily determined with 
the aid of a model or sketch corresponding to 
Fig. 1. After several trials and intensity calcula- 
tions it was found that the best agreement is 
obtained with an arrangement wherein three 
oxygens, marked 0; in Fig. 1, form a ring about 
each gadolinium whereas three other oxygens, 
marked 02 in Fig. 1, form a somewhat closer 
ring around each threefold axis equidistant be- 
tween two gadoliniums. 

Table I shows the results of the final intensity 
calculations leading to a choice of parameters. 
The calculations were made by the usual formula, 


_ 1+cos* 26 
I~ F*. j-—___—_, 
sin? @ cos 6 


In order to be consistent f values for Gd, C, 
and O atoms,’ not ions, were used throughout, 
though it is not to be supposed that the con- 
stituents are wholly in the atomic state. In any 
event, no great change in the calculated intensi- 
ties would result if ionic f values were used 
wholly or in part. 

All visual estimates of intensities given in 
Tables I and II were checked by photometer 
curves. 

For each set of values for the oxygen param- 
eters there are two alternative sets of values for 
the carbon atoms, corresponding to a rotation of 
the formate group through 180° on the line con- 
necting its two oxygens and bringing the whole 
group back into the plane of symmetry. Columns 
3 and 4 of Table I show intensities for a structure 
with three oxygen atoms forming an equatorial 
ring about gadolinium, the carbon lying above 
(column 3) or below (column 4) the level of the 
oxygens as drawn. The general agreement is good 
in both cases, but calculated intensity of 303 


8R. W. James and G. W. Brindley, Phil. Mag. 12, 81 
(1932). 
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TABLE I. Calculated and observed intensities in a powder 
pattern of gadolinium formate. Cr—Ka radiation. 











1 2 3 4 5 6 7 
Indices Intensities 
Due to 
gado- 
Hexa- Rhombo- linium 
gonal _hedral Calculated Observed* alone 
Zz x z x z x 
C 0.12 0.59 0.45 0.86 0.43 0.85 
oO 2 27 2 4 .w® ai 
Ke 36 56 35 S68 SS 
1120. »=— 101 406 426 418 S 415 
1011 100 a 147 139 M 167 
3030 112 44 43 43 98 
0231 aT "th a9 vd M-S = 196 
2240 202 90 81 83 = 65 
2131 ©2071 120/710 144/225 130/715 bie $5195 
1341 212 87 91 1 97 S 1 70 
4150 = 312 42 45\ . 46 56 
4011317 — - _ "- oh ae 
0112 110 47 44 44 M 26 
3251 302 65 73 72> M-S 45 
2022 200 22 22 16 VW 22 
3360 303 30 35 35 WwW 19 
1232 211 26 32 40 M 37 
0551 223 16 22 19 VW 17 
2461 313 41 34 30 W-M 32 
3142 301 22 W 32 
6060 422 12) 15 
5161 412 30} 53 M 30} 60 
0442 «= 223 11) 15 
5270 314 25 W 29 
2352 312 30 W-M 29 
4371 403 22 VW 30 
5052 411 22 VVW 16 
0003 111 5 —_ 5 
4480 404 23 zt. 
1671-324 4 = = — s 
4262 402 39 W-M 34 
1123 210 27 W 35 
1562 323 29 W 37 














* Scale of observed intensities: VS very strong, S strong, M medium, 
W weak, VW very weak. 


exceeds that of 211, contrary to observation. 
This discrepancy is only slight when the carbons 
are below (column 4). It is removed and good 
general agreement extended to all lines by 
slightly tilting the formate groups from the 
positions first considered so that the ring of 
oxygens closest to gadolinium lies at a slightly 
higher level (column 5). (Compare lower part of 
Fig. 1.) 

In all cases the observed intensity for the first 
two or three lines is considerably lower than 
calculation indicates. This is undoubtedly due to 
the high absorption of gadolinium for x-rays. 
On this account intensity observations with Cu 
radiation were disregarded since the absorption 
of gadolinium for that radiation is much greater 
than for Cr radiation. Mo radiation would be 
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absorbed less, but the intensities of lines of low 
glancing angle cannot be closely estimated on the 
Mo patterns. Since any correction would have to 
be largely arbitrary none was attempted. In all 
other respects the final agreement, columns 5 
and 6, is satisfactory. 

A further check can be obtained from the 
calculated and observed intensities for the zero 
layer line of a [111] axis rotation pattern with 
Cr radiation, Table II. The habit of the crystals, 
uniform slender hexagonal prisms, closely ap- 
proaches the cylindrical shape which is ideal for 
intensity estimation in this case. Though the 
observed intensity of 101 again is lower than 
indicated by calculation the remaining agree- 
ment is most satisfactory. This corroborates the 
parameters chosen from the powder intensity 
calculations, namely : 


Gd in 000 


z x 
3 C in xxs; x2x; 2xx 0.43 0.85 
3 O; in xx2; x2x 3 2xx 0.19 0.81 
3 Oo in xxZ; xBx 3 2xX 0.33 0.58 


Nearly one-half of the total scattering powder 
of the crystals belongs to Gd at low glancing 
angles and at higher angles the scattering powder 
of Gd greatly exceeds that of the formate groups. 
As a consequence the intensities are largely 
determined by the Gd. Column 7, Table I, 
shows the powder intensities due to Gd alone. 
It is seen that they show a rough correspondence 
to the observed intensities, but good agreement 
is obtained only upon.taking into account the 
formate groups. Nevertheless, small shifts in 
parameters of the formate group have but a 


TABLE II. Calculated and observed intensities in the zero 
layer line of a [111] axis rotation pattern of gadolinium 
formate. Cr— Ka radiation. 








Intensity 





Indices Intensity due to 
Hexagonal Rhombohedral Observed* Calculated Gd alone 
1120 101 S 45.3 45.2 
3030 112 W 8.3 19.0 
2240 202 VS 18.5 14.3 
4150 312 M-S 13.5 16.2 
3360 303 M 11.6 6.5 
6060 422 W 4.5 | 
5270 314 M 10.0 By 
4480 404 M 9.8 a 
7180 325 VS 31.4 26.0 








* Conventions for observed intensities as in Table I. 
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slight effect on calculated intensities due to the 
predominant influence of the heavy Gd so that 
no great accuracy can be claimed for the chosen 
parameters. 


DISCUSSION OF THE STRUCTURE 


The dimensions of the formate group corre- 
sponding to the chosen parameters may be com- 
pared to those found by electron diffraction in 
gaseous formic acid® and in crystalline sodium 
formate’ as follows: 


Formic Sodium Gadolinium 
acid formate formate 
Bond angle 125° 124° Zt” 
C—O, 1.27A 
Cc-O 1.29A 1.27A , 
C—O: 1.33A 


In view of the different position of O,; and Oz 
in Gd formate it seems improbable that there 
can here be “complete resonance’’ of the two 
C—O bonds as previously found in other struc- 
tures. 

Other interatomic distances are: 


O; to Oz 2.26A O; to O; 3.96A O; to Gd 2.30A 
in same group in ring 
O; to Oz (upward) 2.67A Oz to Oz 2.63A Oz to Gd 2.50A 
in next group in ring 


(downward) 3.01A 


From this it is possible to make an estimate of 
the radius of the Gd ion, since it has been shown 
by crystal structure determinations of numerous 
other compounds that the nearest approach of 
carboxyl oxygens to adjacent metal ions closely 
approximates the sum of the ionic radii. Assign- 


® Pauling and Brockway, reference 6. 
10 Zachariasen, reference 6. 
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2.634 


398A 











Fic. 2. Clinographic projection to show the ninefold 
coordination of carboxyl oxygens about gadolinium in 
gadolinium formate. 


ing to carboxyl oxygen the radius 1.32A, this 
leaves a radius of 0.98A for Gd in Gd formate. 
This value is close to previous estimates of the 
Gd ionic radius." 

The ninefold coordination of oxygen about 
gadolinium is shown in Fig. 2. This is nearly the 
same coordination that has been proposed: by 
Ketelaar™ for Gd and other rare earths in the 
ethyl sulphates, R(C2H;SO,)3-9H:20O. In that case 
nine water molecules are coordinated about the 
rare earth atom. 


11, Pauling and M. D. Shappell, Zeits f. Krist. 75, 128- 


142 (1930). 
2]. A. A. Ketelaar, Physica 4, 619-630 (1937). 
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The Electrostatic Influence of Substituents on Reaction Rates. IV. 


DonNA Price? AND F. H. WESTHEIMER 
Department of Chemistry, University of Chicago, Chicago, Illinois 
(Received January 6, 1943) 


The rates of saponification of ethyl benzoate, ethyl toluate, ethyl anisate, ethyl p-chloro- 
benzoate, and ethyl p-nitrobenzoate have been measured in several alcohol-water mixtures. 
The relative rates of saponification of the substituted benzoic esters have been compared 
with those computed by the Kirkwood-Westheimer theory. 





N a recent investigation,’ the saponification 

rates of the esters of some aliphatic dibasic 
acids were measured. It was found that the 
effect of solvents on the relative saponification 
rates was roughly in accord with the predictions 
of electrostatic theory.4 The present paper deals 
with the effect of solvent on the relative saponifi- 
cation rates of some aromatic esters. In these 
cases, the theory is only qualitatively correct, 
and the work further limits the application of 
the theory to the aromatic series.® 


EXPERIMENTAL 
Materials 


The esters were prepared in the usual manner. 
The liquids were purified by vacuum distillation ; 
the p-nitro derivative by recrystallization from 
70 percent alcohol. The constants were: 


b.p. 99°/10 mm, 

Ethyl toluate b.p. 115°/9 mm, 

Ethyl anisate b.p. 144-145°/10 mm, 
Ethyl p-chlorobenzoate b.p. 123.2-123.9°/10 mm, 
Ethyl p-nitrobenzoate m.p. 57°. 


Ethyl benzoate 


Saponification equivalents for the esters were, 
in every case, within 0.2 percent of the theo- 
retical value. 

The absolute alcohol was dried with mag- 
nesium; the alcohol-water mixtures were pre- 
pared by measuring the calculated volume of 
water into a volumetric flask and filling to the 


1 Samuel S. Fels Fellow in 1942. 

2 Present address: Hercules Powder Company, Wilming- 
ton, Delaware. 

3F. H. Westheimer, W. A. Jones, and R. Lad, J. Chem. 
Phys. 10, 478 (1942). 

4 J. G. Kirkwood and F. H. Westheimer, J. Chem. Phys. 
6, 506, 513 (1938); 7, 437 (1939). 

5 See F. H. Westheimer, J. Am. Chem. Soc. 61, 1977 
(1939) ; 62, 1892 (1940). 
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mark with alcohol. The alcoholic alkali solutions 
were prepared by adding stick potassium hy- 
droxide, which had been thoroughly washed 
with water and then with alcohol mixture, to 
the appropriate solvent. The solutions were 
standardized against potassium acid phthalate 
every three or four days; there was little or no 
change in the titer. The alcoholic alkali remained 
satisfactorily colorless when kept in a glass- 
stoppered bottle. 


Method 


The saponification rates were determined by 
titrating from time to time the amount of 
unreacted alkali present in the reaction mixture. 
The experimental method was the same as that 
used in a previous investigation.* The thermostat 
temperature was 50°+0.02. A measured volume 
of solvent was allowed to come to temperature, 
and a weighed sample of ester, in an open capsule, 
was added to the solvent. The reaction was 
started by adding the potassium hydroxide to 
this solution. The potassium hydroxide, of the 
desired strength, was brought to 50° in a stop- 
pered weighing bottle. The bottle was then 
dropped into the reaction flask and the calibrated 
sampling pipet inserted. Except for the saponi- 
fication of the nitro ester, the sampling pipet 
was filled and allowed to stand for thirty seconds 
with rinse solution in order to minimize the error 
due to change in volume with temperature while 
the sample was being withdrawn. In the case 
of the nitro derivative, each sample served as a 
rinsing for the following sample. 

The solvent and the potassium hydroxide 
solution were measured at room temperature by 
means of a calibrated pipet. The concentrations 
of ester and potassium hydroxide in the reaction 
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TABLE I. Saponification rates at 50°. 


RATES 

















Compound Conc. KOH Conc. ester Conc. KCl “u k in min.~! (m/1)7! 
In 30 percent alcohol 
Ethyl benzoate 0.01015 0.01012 — 0.0102 1.64 
.01015 .01016 — .0102 1.62 
.01015 .01016 0.00757 .0177 1.64 
Average 1.63 
Ethyl toluate .01015 .01014 — .0102 0.740 
.01015 .01015 — .0102 0.775 
.01015 .01016 .00924 .0194 0.752 
Average 0.756 
Ethyl anisate 01015 01015 -- .0102 0.499 
.01015 .01015 — .0102 0.483 
.01015 .01017 01017 .0203 0.493 
Average 0.492 
Ethyl p-chlorobenzoate .004149 .004161 — 0041 4.57 
.004149 .004146 — .0041 4.52 
.004149 .004152 01148 .0156 4.30 
Average 4.46 
Ethyl p-nitrobenzoate .003327 .003322 _— .0033 57.7 
.004149 .004142 — .0041 66.5 
.004149 .004142 .004114 .0083 63.2 
Average 62.5 
In 60 percent alcohol 
Ethyl benzoate .01007 .01032 — .0101 0.641 
.01030 .01032 — .0103 0.657 
.01530 .01523 a .0153 0.636 
.02019 .02019 — .0202 0.641 
.01031 .01031 05147 .0618 0.684 
Average 0.652 
Ethyl toluate .00980 .01005 — .0098 0.326 
.01030 .01032 — .0103 0.312 
.01530 .01524 — .0153 0.299 
.02019 .02019 — .0202 0.301 
.01530 .01530 .0470 .0623 0.304 
Average 0.308 
Ethyl anisate .00980 .01005 — .0098 0.175 
.01030 .01032 _- .0103 0.186 
.01530 .01530 _- 0153 0.184 
.02019 .02019 — .0202 0.179 
.01085 .01086 .0530 .0639 0.180 
Average 0.181 
Ethyl p-chlorobenzoate .00980 .01006 —_ .0098 2.20 
.01030 .01029 — .0103 2.28 
.01530 .01523 _ .0153 2.27 
.02019 .02019 oo .0202 2.29 
.01031 .01031 .01012 .0204 2.52 
.01033 .01033 .01005 .0203 2.19 
.01031 .01031 .03021 .0405 2.38 
.01031 .01030 .04982 .0601 2.38 
Average 2.31 
Ethyl p-nitrobenzoate .00980 .01005 — .0098 40.7 
.01030 .01030 _— .0103 44.8 
.01523 .01524 _ .0152 43.9 
.02019 .02021 — .0202 40.3 
.01031 .01031 .04699 .0573 43.4 
.01031 .01031 .05058 .0609 41.6 
Average 43.0 
In 90 percent alcohol 
Ethyl benzoate .01083 .01082 —- .0108 0.270 
.01072 .01071 — .0107 0.270 
.02091 .02091 .0209 0.265 
Average 0.268 
Ethyl anisate .06051 .06077 — .0605 0.0729 
.06051 .06051 — .0605 0.0716 
.03099 .03102 —_ .0310 0.0715 
Average 0.0717 
Ethyl p-nitrobenzoate .01083 .01086 — .0108 20.6 
.01072 .01072 — .0107 20.6 
.02091 .02091 — .0209 21.2 


Average 20.8 
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TABLE II. Relative saponification rates of substituted esters. 
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Alogk 
Compound 30% alcohol | 60% alcohol | 90% alcohol 
Cale. Cale. | Cale. Cale. | Cale. Cale. 
Obs. R=5.2 R=5.9 | Obs. R=5.2 R=5.9 | Obs. R=S.2 R=5.9 
Ethyl p-nitrobenzoate | 1.58 1.55 0.85 | 1.82 1.60 0.90 | 1.89 1.70 1.00 
Ethyl p-chlorobenzoate | 0.44 0.65 0.35 | 0.55 0.65 0.35 | 0.644 0.70° 0.40° 
Ethyl toluate —0.33 —0.15 —0.05 | —0.33 —0.15 —0.10 | —0.33% —0.15° —0.10° 
Ethyl anisate —0.52 | —0.56 | —0.57 











@ Kindler’s data for 87.83 percent alcohol at 35°. _ 
6 Computed for 87.83 percent alcohol at 35°; see reference (5). 


flask were calculated, however, on the basis of 
the volume at 50°. 

The authors wish to thank Mr. L. Frazier for 
his assistance with the analyses. 


RESULTS 


The rate constants were determined graphi- 
cally and were accurate to about 3 percent. 
Several experiments were tried with varying salt 
concentrations; as might be expected in the case 
of a reaction involving uncharged molecules, the 
salt effects were small. The data are summarized 
in Table I. 

For the esters of monobasic aromatic acids, A log 
k has been defined as the logarithm of the ratio of 
the saponification constant of the substituted ester 
to that of ethyl benzoate at the same tempera- 
ture and in the same solvent. The values of 
A log k have been assembled in Table II. The 
relative saponification rates at 50° in 90 percent 
alcohol are comparable to those obtained by 
Kindler® in 87.83 percent alcohol at 35°. The 
calculated values of A log k, given to the nearest 
0.05 logarithmic units, have been obtained using 
the equations and tables previously published.* ‘ 
The values of the dielectric constant of the 
alcohol-water mixtures were found by extrapo- 
lating Wyman’s’ data from 40° to 50°. The 
values used were 24.0, 37.1, and 53.8 for 90, 60, 
and 30 percent alcohol, respectively. 

The computation involves as a parameter the 
value of R, the distance between the dipole of a 
substituted benzoic ester, and the negative 
charge contributed to the activated complex by 
the hydroxide ion. The computations were made 


6K, Kindler, Ann. 450, 1 (1926). 
7]. Wyman, J. Am. Chem. Soc. 53, 3292 (1931). 


using two different values of R. The first value, 
5.2A, is the one which most closely reproduces 
the data obtained by Kindler®* for a large 
number of esters. When this value of R is used, 
the theory roughly accounts for the solvent 
effect in the saponification of the nitro and 
chloro ester. The observed solvent effect is 
somewhat greater than that calculated; this 
result had been foreshadowed by the work of 
Kilpatrick and Mears® on the effect of solvent 
on the ionization constants of substituted benzoic 
acids. On the other hand, the value of 5.9A for 
R, is clearly inadequate, and this is the value 
which best fits the data for the ionization 
constants of the substituted benzoic acids. It 
had previously been noted’ that in the aromatic 
as opposed to the aliphatic series, the values of 
A log k for saponification are considerably larger 
than the values of Alog K for ionization. No 
explanation for this fact is at present available. 

It should be noted further that the change in 
solvent does not have a large effect on the 
Alog k for the saponification of ethyl anisate. 
This fact is quite reasonable. The slow rate of 
saponification of the ethyl ester of dimethyl- 
aminobenzoic acid has been shown to be due to 
a resonance effect,® and it is to be expected that 
the same would be true for ethyl anisate. The 
stabilization of ethyl anisate can result from 
resonance involving an ionic structure as well as 
the usual one; the change of solvent from one 
rich in water to one rich in alcohol would tend 
to make the ionic structure less stable with 


8M. Kilpatrick and W. H. Mears, J. Am. Chem. Soc. 
62, 3047 (1940). 

9F. H. Westheimer and R. P. Metcalf, J. Am. Chem. 
Soc. 63, 1339 (1941). 
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respect to the neutral molecule. This change 
could only decrease the resonance energy. While 
this effect would probably be very small, it 
would be anticipated that the resonance effect 
would be magnified in solvents of high dielectric 
constant, whereas the electrostatic effect is 
minimized 

The constancy of the relative saponification 
rate of ethyl toluate, as well as the comparatively 
large effect of the methyl group, again suggests 


the type of conjugation between methyl and 
phenyl groups postulated by Baker and Nathan.” 
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I 


ie the Journal of Chemical Physics for May, 
1942, there was published a paper by T. 
Alfrey and H. Mark' on a problem in the statis- 
tical mechanics of binary mixtures which had 
been treated previously by R. H. Fowler? and 
his collaborators by an indirect method. The 
aforementioned authors have attempted a direct 
derivation of Fowler’s results; their attempt, how- 
ever, seems to me to lack mathematical con- 
sistency. 

The problem is as follows: Given a lattice of 
N sites, each being surrounded by z nearest 
neighbors. N4 molecules of type A and Ng 
molecules of type B are placed on these sites. 
There result three different types of pairs of 
nearest neighbors, AA, BB, and AB. Let Naa, 
Nes, Nap be their respective total numbers. 
Obviously, 2NaatNap=2Na and 2Nest+Naps 


. =2Nez so that the three numbers reduce to a 


single independent quantity, (say) Nag. Now, 
the configurational energy depends chiefly on 
these numbers. For the calculation of the con- 
figurational partition function one must know in 
how many distinguishable arrangements of-the 


(1942) Alfrey and H. Mark, J. Chem. Physics 10, 303-4 
42). 

*R. H. Fowler and E. A. Guggenheim, Statistical Ther- 
modynamics (Cambridge, 1939), pp. 244-253 (general); 
pp. 358-366 (liquids) ; pp. 576-581 (alloys). 





molecules over the sites the same value of the 
configurational energy can be found or, what 
amounts to the same thing, the same value of 
Naz. The number of these arrangements is 
called by Fowler and Guggenheim? a ‘‘combina- 
tory factor’? and denoted by g(Naz). If g were 
known, the configurational partition function 
would be given by: 


feontig = Zoltan g(Naz) exp (—Nagw/kT). (1) 


Here, Nagw is the value of the configurational 
energy referred to a special zero. Since the direct 
determination of g is difficult, Fowler and Gug- 
genheim treated (1) in several cases indirectly 
(by means of the “grand partition function”’). 
The value of g is given by them merely as a 
“guess” and can be written as follows: 


g~(NaatNan+Nan)! / 


Nap 3 
Naat( ; ') Nes!. (2) 





Here, a factor independent of Nag remains 
arbitrary. (Note that the numerator is inde- 
pendent of Nz.) 

The form of (2) suggests that g is the solution 
of a simple combinatory problem where a given 
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number of elements is to be distributed over four 
groups of Naa, Nas/2, Nas/2, Nes members, 
respectively. The pairs cannot be treated as 
elements, however; they are not independent 
entities, owing to the overlapping of adjacent 
pairs. Fowler and Guggenheim knew this per- 
fectly well. They wrote: ‘This formula .. . 
cannot be accurate because we cannot treat the 
pairs as independent entities. Obviously, as soon 
as we fix the nature of a certain pair of sites this 
limits the possible manners of occupation of 
other pairs which share a site with the first 
pair.” 

Alfrey and Mark propose to prove (2), and to 
avoid the difficulty just mentioned. They group 
the N sites in N/2 pairs of nearest neighbors so 
that each site belongs to one and only one pair. 
This they call a ‘“‘reference set of pairs.’”’ It should 
be noted that this is only a fraction of the total 
numbers of pairs, Nz/2, which may be formed 
between the sites. If the lattice is filled with 
molecules, three different types of pairs will 
result among the N/2 selected pairs; their re- 
spective numbers I shall denote by 144, mgz, Nap. 
To find the number of distinguishable configura- 
tions, these authors apply the rules of combina- 
torics to the selected set of pairs, treating them 
as independent elements. Accordingly, they find 
a number 


N 
= (F1/mastnan'nnn! Java (3) 


Is this really the total number of distinguishable 
configurations of the molecules over the sites 
which belong to a fixed value of the total number 
of pairs, viz., Naz? Here, the authors adopt a 
hypothesis, namely : 


Nap=2nap. (4) 


That is to say, the total number of AB pairs 
should always be proportional to the number of 
AB pairs found in the selected set. The authors 
do not justify (4); they merely say: ‘“This will 
be true in the overwhelming majority of con- 
figurations.” 

I fail to see that this hypothesis is true. It is 
impossible to deduce from the selected set of pairs 
anything about the total set of pairs. The selec- 
tion is arbitrary. In one case it may have been 
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made so as to include more pairs of the type AB 
than in another; the ratio of the number of the 
selected pairs AB to the total number of such 
pairs will therefore differ from case to case. Nor 
is it possible to assume that the hypothesis is true 
for the probable average of the selections since 
ordinary statistics are not applicable here be- 
cause the pairs overlap and are therefore in- 
terdependent. 

If the hypothesis (4) were correct, it would be 
possible to prove the formula (2) quite generally : 
Note that from (4) it follows that N44=2n,,4, 
Nes=2nee and if Stirling’s theorem is applied 
to (3), this equation reduces to the form (2). 
However, the hypothesis (4) is incorrect, and 
therefore so is the proof given by these authors. 

The rest of Alfrey and Mark’s paper repro- 
duces an argument of Fowler and his collabo- 
rators concerning the evaluation of the partition 
function (1) which is incomplete and requires 
additional justification (see Part II), although 
this concerns a point which is usually taken for 
granted. 


II 


In a paper* presented to this journal on Oc- 
tober 14, 1941, I treated the same subject 
(among other items). I pointed out a special 
type of mixture for which it is possible to prove 
Fowler and Guggenheim’s formula (2). As far as 
I can see, this is the only case where this can be 
done. In view of what has gone before it will be 
useful to repeat my argument. 

The type of mixture adopted whose physical 
meaning will be discussed later has the following 
structure: Let the molecules be gathered in 
clusters and let each cluster be widely separated 
from the others. Let each cluster consist of a 
central molecule and a single layer of surrounding 
outer molecules. The latter’s number, 20, is 
taken to be so small that the mutual distances 
of the outer molecules exceed their distances from 
the center. Under these circumstances the only 
pairs of closest neighbors are those formed be- 
tween the centers and the outer molecules. It is 
obvious that these pairs can be treated as inde- 
pendent entities and hence a proof of (2) be- 
comes possible. 


3 F, Kottler, unpublished paper. 
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The proof is as follows: Let us consider the 
sites of the molecules and fill them with molecules 
A and B. This must be done in such a manner 
that no ordered structure results, that is to say, 
there should be no preference of either type of 
molecules for central or outer sites. 

We begin with filling the central sites. There 
shall be, say, Na’ molecules A and N3z’ molecules 
B on them. Consequently, there remain 29N 4’ 
molecules A and zoN,’ molecules B for the outer 
sites. We may treat separately those clusters 
with A’s at their centers and those with B’s at 
their centers. There are 2)N,4’ outer sites in the 
first clusters altogether. Let X; be the number of 
B’s placed among them; the remaining number 
zoN4’—X,, must, therefore, be filled with A’s. 
The number of distinguishable ways in which 
that can be done is: 


ga= (zo 4’) WX, '(Zol a’ —X1) t. 


In a similar manner we obtain for the number of 
different ways to fill the outer sites of the clusters 
with B’s at their centers: 


gp= (Z0Np’)!/X2!(zoNe’ —X2)!, 


where X,2 is the number of A molecules placed on 
these outer sites. The total number of different 
ways to fill the assumed structure with a mixture 
AB is therefore, 


g=Zgate* (oN 4) '(2oVp')!/ 
X1!(zoN a’ —X1)!Xe!(zo Ng’ —X2)!. (5) 


Now, we must have X,= Xo, since there cannot 
be preference as to the kind of site for any 
member of a pair AB. Obviously, 


X1=X2= Nazp/2, 20Na’—Xi1= Naa, 
zoN pn’ —X2= Nps, 


and hence (5) reduces to the form of (2). 
Having thus proved Fowler and Guggenheim’s 
formula for this special mixture it is easy to 
supply the required justification for their evalua- 
tion of the partition function (1) briefly men- 
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tioned in Part I. They simply replace the sum 
in (1) by the greatest term to be found in this 
sum and remark merely: “. . . the error can as 
usual be shown to be negligible.”’ The justifica- 
tion is as follows: With the form (5) for g, the 
series (1) becomes a finite hypergeometric series. 
The sum can therefore be expressed by a complex 
integral.* This integral, in turn, can be evaluated 
asymptotically by the method of the saddle 
point.5 The result so obtained coincides with the 
greatest term in series (1). 

It remains to discuss briefly the physical 
meaning of the structure adopted for the mix- 
ture. It has certainly nothing to do‘ with a 
crystalline lattice. Therefore, the above reasoning 
cannot be applied to the crystalline state.* But 
it seems to me that a limiting stage of the liquid 
state may exist where it applies. As a matter of 
fact, a dense gas in the transition to liquefaction 
may consist of the clusters as described. If the 
density of the liquid increases, these clusters 
will acquire more than one outer layer and their 
distances will decrease. Moreover, the coordina- 
tion number will increase from the low value 2» 
to the value for closest packing. There is no 
experimental evidence, so far as I can see, 
against the assumption of such a limiting stage. 
For liquid mixtures, formula (2) may, therefore, 
be regarded as a useful approximation which is 
accurate for the assumed limiting stage of the 
liquid. Unfortunately, our knowledge of the 
liquid mixtures is not sufficiently advanced to 
permit an examination of the refined theory 
based on formula (2). 


4E. T. Whittaker and G. N. Watson, A Course of 
Modern Analysis (Cambridge University Press, Cambridge, 
1940), p. 292. 

5R. Courant and D. Hilbert, Methoden der Mathe- 
matischen Physik 1 (J. Springer, Berlin, 1924), pp. 430 
et seq. 

*Formula (2) for a crystalline mixture appears to be 
incapable of direct proof and is probably not true. Some 
indication of this may be deduced from the numerical 
failure of the deductions based on the formula (2) or its 
analogue, in the explanation of certain anomalies of the 
specific heats of alloys (reference 2). 
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The Electrical Capacity of the Double Layer 
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N a recent article,! J. J. Bikerman has extended the 
theory of the diffuse double layer to high potentials and 
high concentrations by taking into account, (1) the ‘‘satu- 
ration”’ effect which relates the maximum possible concen- 
tration of ions in the double layer to the size of the ions 
and the space which is available, and (2) the electrical 
energy resulting from the polarization of hydrated ions and 
polar solvent molecules. It becomes desirable, therefore, 
to define the range of “low” potentials and ‘‘low’’ concen- 
trations in which the equations of a recent paper? of the 
present authors are reasonably valid. This can be done 
only in a very approximate fashion since the evaluation of 
the factors (1) and (2) requires a knowledge of the volume 
of the hydrated ion and its effective moment. However, at 
least the order of magnitude of this range can be esti- 
mated.* 

(1) The saturation effect should not be the limiting 
factor as long as the total ion concentration is less than 
12 moles per liter. This is less than the ion concentration 
in a saturated solution of sodium chloride and corresponds 
to the condition, 2C® cosh 2¥=12, which gives Y=60 milli- 
volts if C°is 1 mole per liter and y=290 mv at 1074 molar. 

(2) The polarization effect is more difficult to estimate 
since it involves the difference between the polarization of 
a hydrated ion and that of an equal volume of solvent 
molecules. However, the ion should displace approximately 
one water molecule and exert an undetermined influence 
on a few others distributed spherically with respect to the 
ion. This effect should not be important if the electrical 
energy of an ion (ey) is much greater than the electrical 
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energy of a water dipole (m(dy/dx)) where m is the 
average moment of the dipole in the direction of the field. 
For a water molecule (assumed rigid, with a moment of 
1.84 10-8 and with free rotation) the value of the ratio 
r (=m(dy/dx)/ep) in 1 molar solution is 0.014 at y=50 mv 
and 0.064 at Y=100 mv. In 10-4 molar solution r=3 X 10-5 
at Y=150 mv and r=3X10-3 at y=260 mv. 

Therefore, the neglect of these factors should not greatly 
reduce the accuracy of the equations derived in the previous 
paper? in the experimentally determined range of ¢-poten- 
tial values. ¢-potentials greater than 100 mv have rarely 
been observed and such values have been found only in 
very dilute solutions where ¢ and y are not very different 
(Eq. (3)). 

Bikerman’s equation. for low potentials, on the other 
hand, involves the assumption that sinh 2)=2y, and, 
therefore, (allowing a 10 percent error) use of this equation 
should be limited to the range of potentials where sinh 2y 
=1.10(2y). This corresponds to the condition y=20 mv 
at all concentrations. At such low potentials (1) and (2) 
are negligible and, therefore, his equation for the capacity 
of the double layer is the classical one, C’ = Dx/47, values 
of which are given in Table II.? This gives a ‘‘thickness”’ 
of the ion atmosphere (1/x) which is independent of the 
electrical condition of the surface. Actually this thickness 
decreases rapidly with increasing values of the surface 
potential (Eq. (12)). 

Bikerman’s equation for high potentials may be applied 
to the conditions existing near the surface where higher 
field strengths lead to “‘saturation” of the surface volume 
and deformation of ions and molecules. It should not be 
extended to low potentials (or to greater distances from the 
surface) since it assumes that one of the ion species is 
excluded from the outer double layer. From the standpoint 
of the present equations? for “‘low’’ potentials its im- 
portance lies in the region between the surface of the solid 
and the y’ plane. For ‘‘high”’ potentials this may be greater 
than a single molecular thickness. It has long been recog- 
nized that such a region is frequently required to account 
for a large difference between the total potential difference 
across a phase boundary and the ¢-potential. For example, 
if the ¢-potential of mercury were equal to the total poten- 
tial difference between mercury and a normal solution, 
a mercury droplet detached from a calomel electrode should 
be expected to have a cataphoretic mobility corresponding 
to a ¢-potential of over 500 mv. Such high cataphoretic 
mobilities have not been observed even in dilute solutions. 

1J. J. Bikerman, Phil. Mag. 33, 384 (1942). 

2 W. G. Eversole, and C. R. Estee, J. Chem. Phys. 11, 63 (1943). 


* The symbols and equation numbers of this note are those of 
reference 2. 








